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NOTE ON A PROBLEM IN ELECTROSTATICS . 


By E. P. ADAMS (Princeton) 
[Received 11 February 1939] 


1. THE ingenious method used by Richmond (1) for getting the 
solution of problems in electrostatics connected with a grating of 
rounded bars may be applied to the problem of such a grating placed 
midway between two infinite planes at the same potential. 

We begin with a grating of cylinders of rectangular cross-section 
as shown in the figure which represents the z-plane transverse to the 





















































axes of the cylinders. The distance between the two parallel planes 
is 2h, and 2d is the distance between the axes of consecutive cylinders, 
each of dimensions 2a along the z-axis and 2b along the y-axis. On 
account of symmetry we need consider only the area within the 
dotted lines. This area is transformed to the upper half of a ¢-plane 
by the Schwarz-Christoffel method, the real values of ¢ at the corners 
being O0< 1<s<n<m< oo. The differential equation is 

dz A(ji—s}t 

dt  {t(t—1)(t—n)(t—m)}* 
Following the method of Richmond, we put 

A(t—s)t = P(t—n)'+ Q(t—1)!, 
and get 
dt 


dt 
~<e J fut—1)(t—myji* @ J e—ny(t—myi* > 


where B is a constant. In the first integral put 
dn?w 


with a modulus k determined by 


= 1—-, 
m 


3695.10 R 








242 E. P. ADAMS 
and with corresponding values of ¢ and w given by 
t= 0 1 m 6) 
w=tK’ 0 KK K-+K’. 
For t = n, since 1 < n < m, w = Up, where uy is real and lies between 
0 and K, and is determined by 
1 
dn*u, = -. 
n 


In the second integral put 
n 


= — 
dn*w, 
with a modulus k, determined by 


9 n 
=1—%. 
m 


Corresponding values of ¢ and w, are given by 
t= 0 n m 00 
oe r r 4 sere 
w,=tK, O KK, K,++#K;. 
For t = 1, w = ivy where vy is real, lies between 0 and Kj, and is 
given by ke 
2/9’) - 1 
dn?(v>, k;) = ey 
We can now write 
, 1 — 
z= £f w+ Qu 17 B, 


and the constants are determined from the table, 
[t= 9 LL n- om o.0) 
<= th ha dad d+th 
w=1K'’ 0 u K K-+iK’ 
w,=—tK, i O K, K,+#K, 
On substituting these values we get 
B= @, d= P’'K+QK,, b= Qn, 
a= Pm, h = P'K'+Q’K;j. 
If we take a = b so as to have an approximation to circles, we have 


p= 


+ (K+ “*K,), 


Ug Vo 


+(x’+“2K;), 


Uo Vo 








ON A PROBLEM IN ELECTROSTATICS 

and Up, UV) are determined by 
k’2 k2—k2 

: sna, k) = —_+, 
ki? ( 0 ) kk;? 
ki k?—k? 
k2’ k.3 © 
By assuming values for k and k, (< k) the ratios d/a and h/a may 
thus be found. 

On the cylinder, 1 < ¢ < n, so that w = u with 0 < u < K, and 
w, = ww with O << v < wv. Put 

t = 1+(n—1)A, 

so that A= 0, 1 for ¢=1, nm. Then the equation of the positive 
quadrant of the cross-section of the cylinder may be written 


dn*(up, k) = 


dn*(v9, ky) = sn*(vp, ki) = 


x 


me 


> 


u 
ec % &© & 


and wu, v are given by 


2 k’ 
dn (u, k) = k2—(k2— Ky’ 


; je2 
dn*(v, kj) = Be R 
For a circle we should have 
OF cet 
a® a? ; 
Milne-Thomson’s tables (2) are useful for rapidly getting a numeri- 
cal solution. As an example, take 
if == §'* = 05; H=01, £* = 04. 
Then we find Uy = 1-3821, Vp = 1-6327, 
and d/a = 2-3291, h/a = 2-9206. 
For the shape of the cylinders we find 
A=@ 0-1 0-3 0°5 0-7 0-9 
u=0 0:4003 0-6969 0-9098 1-0968 1-2819 
v = 16327 1-5677 11-4239 1-2505 1:0196 0-6232 
z/ja =0 0:2896 0-5043 0-6583 0-7937 0-9276 
yja=1 0-9603 0-8722 0-7659 0-6245 0-3818 
rja=1 1-0030 1-0072 1-0097 1-0095 1-0031 
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Thus for all practical purposes these cylinders can be regarded as 
circular. 

If k =k, and k’ = k,, with k, < k, we get h/a = d/a, uy = %, 
u =v. The cylinders are now symmetrical not only with respect to 
the x- and y-axes, but also with respect to axes inclined to these at 
an angle of 45°. Thus with k? = 0-9 we find 

h/a = d/a = 1-87738, 

and the greatest value of r/a, for A = 0-5, is r/a = 1-034. For 
2 — 0-8 we find hia = dja = 2-2598, 
and the greatest value of r/a is 1-015. 

For k and k, small, less than about 10-*, we can put 


KX = K, = is, 
ry rs 4 
K’ = log- K,; = log—, 
ky 


ki 

x7 2. — gin2, eine 

sn*uy = sin*ty = 1-73, 
sn?(v9, kj) = tanh*v, = sin?2Wp, 

so that SCC Uy = cosh vp, 


d — 1(Uy+%p) 


a 2Uy Vo 

1 4 u + 
— = —(log—+— log —). 

z( ure 4 
For k = k, = 0, h/a = ow and we have the case considered by Rich- 
mond of a grating of rounded bars in free space. 

2. Let Q be the charge on unit length of each cylinder at potential 
V, the two infinite planes being at potential 0. The complex poten- 
tial is ; 

’ x= b+ip, 
with ¢ as the real potential. The line of force from t = 0 tot=1 
is taken as ¢ = 0. Then the line of force from t= n to t= & is 
ws = 7Q. The equipotential from ¢ = 1 tot = nis ¢ = V, and from 
t=oo tot=0 it is 6 =0. The y-plane is now the interior of a 
rectangle, and to transform to the upper half of the t-plane we have 


dy _ C 


dt {t(t—1)(t—n)}*" 
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1 
dn2w,’ 
with a modulus k, determined by 
2 = m—1 _ ke— ky 

, n kb? 
We then get x = C'w,+D, 
and the constants are determined by 


Put ‘= 


i= © 1 n (ere) 
W,=tK, 0 Ky K,+iky 
x= 0 V V+irQ inQ. 
We then get x= we wyt V,~ 
2 


Q_ , 
and V "ake 
The last equation gives the ratio of the charge on unit length of 
each cylinder to the difference of potential between the cylinders 
and the planes. With 
_7K _7Ks 
Z=¢€™, G=e™®, 


Q 1 1 


ee =— —] 
Vs ilogi/qg._ +7 - 


we get 


1 

“ 

In the example considered, with k? = 0-5, k} = 0-1, we get 
kz = 4/9, 


and ¥ = 0-3025. 

By giving different values to the ratio b/a the exact solution may 
thus be obtained for a grating of cylinders of a variety of shapes 
placed midway between two parallel planes at the same potential. 

3. If we solve the problem of an infinite series of line charges Q 
at a distance 2d apart, placed midway between two infinite planes 
at the same potential at a distance 2h apart, and assume that 
cylinders of radius a with the line charges as axes are the equi- 
potential surfaces V, we shall find 


Q d 
vo i/ (210e aKa) 
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with a modulus determined by 


For the example considered, with h/a = 2-9206, d/a = 2-3291, we 
find 
Q = 0°2995, 


V 


which is only about 1 per cent. less than the value we found for the 


approximately circular cylinders (0-3025). 

4. The numerical results obtained by Richmond may be used to 
solve the problem of a single cylinder midway between two infinite 
planes at the same potential. This is the condenser problem that is 
solved for a circular cylinder by Love (3) in a recent paper in this 
Journal. It is readily found that the capacity S is given by 

K 
S = >? 
aK 


where the modulus is k = sin, 
and cosh b = sec, = a ; 
z(b+) 
With the values of 8 given by Richmond (except for the last for 
which I find b = 3-9045, B = 87° 41’ 25”) we get the following table: 
a/h = 0-1 0-2 0-5 0-7 
S = 0-1965 0-2702 0-5433 0-9318 
C = 0:1965 0-2701 0-5381 0-8684 
S’ = 0°1965 0-2701 0-5349 0-8339 
The third line gives the values of the capacity C, which were found 
by Love for accurately circular cylinders. The last line gives the 
values of the capacity if we assume that the circular cylinder of 
radius a with a line charge Q midway between the planes as an axis 
is the equipotential V. The capacity in this case may be shown to be 


SS = i/(2 log} 
77 


REFERENCES 
1. H. W. Richmond, Proc. London Math. Soc. (2), 22 (1924), 389. 
2. L. M. Milne-Thomson, Die elliptischen Funktionen von Jacobi (Berlin, 1931). 
3. A. E. H. Love, Quart. J. of Math. (Oxford), 9 (1938), 246. 








TRANSFORMATIONS OF ENUMERABLE SETS 
WHICH ARE DENSE IN AN INTERVAL 


R. COOPER (Belfast) 
[Received 28th March 1939] 


1. ONE-ONE bicontinuous transformations of enumerable sets which 
are dense in given intervals have been studied by several writers. 
Among their results are the following: 

(i)* Let #, F be enumerable sets, each dense in Euclidean space 
of nm dimensions. Then there is a one-one transformation P = f(p) 
in which P,p are points of H,F respectively such that f and its 
inverse both satisfy a Lipschitz condition, i.e. there is a constant K 
such that 1 _ p(P,B) 7 

K ~ p(Py P2) 
where (Pi, M2), (P2, 2) are pairs of corresponding points and p(a, b) 
denotes the distance between the points a, b. 

(ii) W. Sierpinski} says that a set F is smaller-(N), or smaller in 
von Neumann’s sense, than a set # if there is a function f(x) defined 
in EF which sets up a one-one correspondence between FH, F and 
satisfies the condition that |f(x)—f(y)| << |w—y| for each pair of 
distinct points x, yof H. He proved that, if Z, F are linear enumerable 
sets contained in and dense in intervals of lengths a,b and if b < a, 
the set F is smaller-(N) than EZ. It is obvious that F cannot be 
smaller-(N) than # if b > a. 

2. In this paper we are concerned primarily with the case not 
covered by Sierpinski’s result, ie. when the intervals are of equal 
length. Our result when combined with a linear transformation 
throws light on the transformations which can be constructed in the 
case of unequal intervals. If either of the intervals is unbounded, 
separate consideration is needed, but it will be sufficient to mention 
briefly the differences from the argument of the primary case. In 
all cases. it appears that a very smooth relation can be established 
between any pair of enumerable sets dense in intervals. The following 
theorem gives the principal result. 


* P. Urysohn, ‘Sur les multiplicités Cantoriennes’, Fund. Math. 7 (1925), 
29-137, esp. 83-9. 

+ W. Sierpinski, ‘Sur une propriété des ensembles dénombrables denses 
dans un intervalle’, Mathematica, 11 (1935), 222-8. 
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THEOREM. Let HE, F be enuwmerable sets contained in and dense in 
the interval (0,1). Let « (<1) bea small positive number and let N be 
any integer. Then there exists a function f(x) defined forO <4 <1 
which sets up a one-one correspondence between E, F and satisfies the 
following conditions: 

(i) f(x) is convex and increasing in (0, 1); 
(ii) f(x parti O<x<landv<QN; 

= l—e — <I+e for 0<2x<l; 

) |fm(ax)| << « pee: y= 2,3,..,N and O0<2<l. 

We construct this function as the limit of a decreasing sequence 
of functions f,,(~) each of which sets up a one-one correspondence 
between n pairs of points of #,F. If N is unity,} it is possible to 
take these approximating functions as polygonal functions,{ but, 
if on such lines we try to satisfy conditions about the second eaies- 
tive, great difficulties arise, partly on account of the discontinuities 
of the derivative of the approximating functions. We therefore use 
as approximations functions which have no discontinuities among 
their derivatives up to the Nth order. 

The construction arranges the sets H, F as sequences {x,}, {y,}, 
and we are able to define a function f,(x) to satisfy the following 
conditions in (0,1): 

(i) f,(#) is convex and increasing; 
(ii) f(x) exists and is = forv < N; 
= Tet 3-*- < f(x) < 1+e—2-"—1¢; 
iv) |f@(x)| << (l—2-"- a 3c va: 
o the x er un of f;(x) is not zero; 
(vi) f,(0) = 0, L(ij=i1, Findj=—y, (r= 1, 8...., 9). 


We also are able to arrange that, for 0 < # < 1, 


5 id 9-—n-2 
x) Fn+al x) )| <2 €, 


and fP@)—fesale)| <2" (1<v <Q), 
From these it follows that the sequences {f,(x)}, {f(x)} (v<N) 
are uniformly convergent in 0 < x < 1 and f(x) = lim/f,(7) has the 
properties stated in the theorem. 


* The derivatives at 0 and 1 are right- and left-hand derivatives. 

+ This was the case I originally considered, and I am indebted to Dr. J. 
Gillis for suggesting that it might be possible to impose condition (iv). 

{ This is on the lines of Urysohn, loc. cit. 
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3. We now give the details of the construction. The curve 
y = f,(x) is not required to pass through any points (z,,y,). We 
therefore choose it to be a circular arc. 
To construct the function f,, ,,(~) from f,,(x) we require the following 
lemma: 
LeMMA. Given any positive number 5, a function ¢(x) can be 
defined which will satisfy the following conditions in an interval («, B): 
i) Ja) = 0, 0< (8) <8; 
(ii) $a) = $B) =0 for 1<v<VQ; 
(iii) O< g(x) <8 for a<a<B; 
(iv) |¢™(a2)|<8 for 2<v<cN,a<2x<B; 
v) d(x) is continuous. 


For this it is enough to observe that, if k is sufficiently small, 


2) =k f (t—a)*(B—t) dt 


satisfies the conditions. 

Since the sets #,F are enumerable, we can arrange them as 
sequences {£,}, {n,} say, The curve y = f,(x) passes through the 
points (x,,y,) (7 = 1, 2,...,2) whose coordinates are certain terms of 
the sequences {€,}, {n,}. ; 

When » is odd, we take z,,,, to be the first term of {€,} which is 
not one of the numbers 2, 29,..., x,. The numbers 2j,..., %, divide 
(0,1) into n+1 intervals; let (a’,2”) be the one which contains 
X43, and let y’ = f(x’), y” = f(x”). We choose for y,, ,, the first term 
in {,} which lies between y’ and y” and is such that there is a 
function f,,,,(7) with the following properties: 


(i) frsi(v) = f,(x) except for x’ < x < x”, and 
fnsa(®) < f(x) for 2 << a< 2"; 

(ii) fnsa(@nsas) = Ynsai 

(iii) f,(z) exists and is continuous for 0 < x < 1; 

(iv) fnr4a(v) is convex and strictly ipereasing; 

be ) FEL (ae) = f(x) at x’, 2,4, and 2” (l<v< N); 

vi) |f) (2) —f@(x)| < 8, (0 <v < N) where 5, is a prescribed 
Be which is less than 2-"-*e. The lemma shows that this is 
possible provided that 5, is less than the lower bounds of f,(z) 
and f, (2). 

When n is even, we choose y,,,, as the first term of {n,} which is 
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not one of y;,..., y,, and then xz, ,, to be the first available term of 
{é} so that the same conditions are satisfied as in the case of n odd. 
This completes the construction; it is obvious that the sequences 
fx,}, {y,} contain every term of the sets EZ, F. 

4. We see at once that, if the given sets ZH, F are contained in and 
dense in different intervals, there is a function transforming E into 
F which has the same properties as the function of the theorem 
save that instead of the condition (iii) its derivative lies in a pre- 
scribed neighbourhood of the ratio of the lengths of the intervals. 

5. Unbounded Sets. If both the sets #, F are dense in (—0o0, 0), 
we can construct a function with the properties of the one in the 
theorem save that they hold for all values of 2 instead of in (0, 1). 
Only slight modifications of the above argument are needed; if 
k is sufficiently small, we can take for a first approximating function 

folx) = Hx+y(1+h*2*)}, 
whose graph is one branch of a hyperbola with asymptotes having 
gradients $(1-L-k). 

We could similarly require the derivative to differ by not more 
than e from any selected positive number instead of from unity. 

If the intervals containing H, F are (a,00), (b,00), we can satisfy 
the same conditions as in the case of the intervals (—00, 00), (—00, 00). 

6. Further differences appear in other cases. First suppose that 
the intervals containing #, F are (—0oo,00) and (a,b), where a,b are 
finite. Then we can select arbitrarily a point c and a small number e, 
and show that there exists a function f(x) such that y = f(x) sets up 
a one-one correspondence between x in EF and y in F and such that 
(i) f(x) is strictly increasing; 

ii) f(x) is convex in (—oo,c) and concave in (c, 00); 
ii) f(x) exists and is continuous for all x for v < N; 


rf 
\ 
(i 
(iv) |fM(z)|<e (laxv<WQ). 
The argument is similar to those above. We can construct a first 
5 
approximation satisfying the conditions (i)—(iv) with 4e in place of e, 
e.g. if k is sufficiently small, we can take 


_ b+a_b 


fo(x) = —5-+ —“tan-k(e—c). 


If c is not in E but $(b+a) is in F, we take y, = }(b+a). If c is in 
E but }(b+-a) is not in F, we take x, = ¢ and reverse the rules above 
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about odd and even suffixes. Ifc is in Z and 4(6+-<a) is in F, we take 
x, =cand y, = 3(b+a) and f,(x) = f(z). 

In addition to using the lemma as above we have to apply it as 
follows. If x, ,, is to the left of all the points x,, %,...,2,, we take a 
number X < 2,,, and apply the lemma to the interval (X,z,,,,), 
because f,,,, and f,, must tend to the same limits as z > +00 and —oo. 
The case of x, ,, to the right of all of z,,...,z,, is treated similarly.* 

7. If F is in an interval (a,00) and E in (—o,00), we can find a 
function which transforms EF to F, satisfies the conditions (i), (iii), 
(iv) for f(x) in § 6, and is convex in (—oo, 00). To see this it is sufficient 
to note that a suitable first approximation is 

a+ka+./(k®x?+ 1) 
if k is small enough. 

8. If Z isin (a,0o) and F in (b,c), there is a function f(x) transform- 
ing £ into F which satisfies 

(i) f(x) is concave and increasing for a < 2; 
(ii) f(x) exists fora < x and v < N; 

(iti) |f(x)| << « fora <2. 

For this we need only note that 
(c—6) 


~ tan-"k(a—a) 
7 


2 
b+ 


is a suitable first approximation if k is small enough. 

9. The only cases not examined above are those in which one of 
the intervals is of the type (—oo,a), but a reflection reduces these 
to cases already dealt with. 

It is possible to establish similar results for sets in space of any 
number of dimensions, but the results and methods are so similar 
to those used above that it is unnecessary to state them. 

Our results show that the possible degree of smoothness of the 
transformation function can give no method of assigning any relation 
of relative density to a pair of enumerable sets. 


* This step was not needed in § 5, because there, in the case of 
Ta <% (¢ = I,.... 2), 


we could take fp,1—f, constant in (—, %p_1). 
+ I am indebted to the referee for some suggestions which have enabled me 
to make some points clearer. 











ON SELF-RECIPROCAL FUNCTIONS 
By B. MOHAN (Benares) 
[Received 1 May 1939] 


1, A Few years back I proved some theorems concerning self- 
reciprocal functions. The object of this note is to give some more 
general theorems of the same type and to derive some new self- 
reciprocal functions therefrom. 

I shall say that a function is R, or —R, according as it is self- 
reciprocal or skew-reciprocal for J, transforms. 

Hardy and Titchmarsh have defined a class of functions A(«, a), 
where 0 < a < 7,a < },as functions (i) which are analytic functions 
of x = re? regular in the angle A defined by r > 0, |0| <a, and 
(ii) which are O(|x|-¢-*) for small x, and O(|a|¢-1+*) for large x, for 
every positive 6 and uniformly in any angle |@| < a—7 < a. 

I take the following theorem from Hardy and Titchmarsh (7) 
concerning this class of functions: 

A necessary and sufficient condition that a function f(x) of A(a,a) 
should be R,, is that it should be of the form 

; ct+io 
f2) = = | 2'T'(}+-dv+ 38)p(s)a-* ds, (1.1) 
27 S ? 

c—ia 
where us(s) is regular, and satisfies the equation 

¥(s) = ¥(1—s) 
in the strip a<o< l—a; 

b(s) = Ofett7—a+mitlt (1.4) 
for every positive » and uniformly in any strip interior to (1.3); and 
c is any value of o in (1.3). 

2. THrorEM I. If 


b+ia 
1 
Pea) = 5 [ 2T+Hn+ $9) b+dy-+Jo)o(6)2-* as, 
b—iaw 


where w(s) satisfies (1.2) in the strip (1.3), and 


w(s) = Ofettr—20 +27) tits 
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and k+io 
fle) =f DEH Hd) +4 + P—do)x(o edo, 
&-to (2.3) 
where x(s) satisfies (1.2) in the strip (1.3), and 
x(s) = Ofeiz—a+nit} (2.4) 


then the function g(x) = [P (xy) f(y) dy 
0 


is R, and belongs to A(3a, a). 
By (2.1) we have 
1 rc) b+io 
ate) = se [ fu) dy | 2T+En + B+ y+ do)a(o) ey) do 
0 b-io (2. 5) 
b+io co) 
aT b+ Ie bOT(E+ b+ isde(e)e* ds | yey) dy. 
b—iaw 0 (2.6) 
To justify the inversion of the order of integration we observe that 
the modulus of the repeated integral is not greater than a constant 
multiple of 


«4 f 2>|P(4-+-du-+ 40+ fit| |P(4+ 4+ 30+ Bit) | |w(s)| dt x 


— 2 


x [ y*Ify)| dy. (2.7) 


The inversion of the order of integration will be justified if the 
two integrals separately exist. 
Now, by virtue of (2.2) and 
IP(B+4it)| ~ Cerne |2, (2.8) 
the integrand in the t-integral is substantially the same as 
etal |¢|—t+4e theta |g |—t+hv +hoedar—2a+2n)iE\ 


i.e. the same as e—%a—itl |g |b tbe thy +0, 


Hence the f-integral exists for a—n > 0. 
Also, by virtue of (2.4) and (2.8), the modulus of the integrand in 
(2.3) is not greater than a constant multiple of 
etal |g] bette — del |g |b tte dell |g |b +bA—theer— at lly 
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which is the same as 
en (a—mylll |g | tte tiv +tAt+ akg —k, 
Therefore, the integral in (2.3) is uniformly convergent for « > 7. 
Also, f(x) ~ |x|-*, and since k may be taken as near b or 1—b as we 
please, it follows that the y-integral in (2.7) also converges. Hence 
the inversion of the order of integration in (2.5) is justifiable. 
Now, by Mellin’s inversion formula (6), we have, from (2.3), 


[ e-fte) de = HE h+ aut do) d+ 49) G+ — Is) x00). 


Thus 
| a-of(ar) de = BVT (E+ du — 4s) (3+ —49) (E+ +48)x(6), 
0 

since x(s) satisfies (1.2). Therefore, from (2.6), we have 


b+ia 
a 


where 
f,(8) = Vi $8)P'($+34—}38) x 
x P'($+ $v+ 4s) 0 (3+ 40—$s)w(s)x(s). 
Now y#,(s) satisfies (1.2), and, by (2.2), (2.4), (2.8), 
ob, (8) cam Ofe-m"Al +(}7r—2.x+2n)|t| (37 -—ax +m) It } 
= Ofett—Sa+3n)t} 


It follows from (1.1) that g(x) belongs to A(3«,a), and is R. 


3. Example I 
I start with the particular case of a formula of Goldstein (5): 
Be ae P(-+-m+4)l(l—m-+4) 
yl—le—ia ‘ ree 2 
a!—le-izpf,, (x) dx P—k-b1) 
where /+m-+4 > 0. This may be thrown into the form 
r 2-1—42" Ed P(/-+-m-+-43)(l—m-+4) 
[« e Wi m( 42") dz = MIP (I—k-L1) , 


0 


> 
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where 1+m+}4>0. Putting 1 = }s+}u+h)—}, m= }(y—»), 
= $+}4+})—)H, we get 


oO 


I as tah te a sty tau 2") Or 
: ptevtuste a ECE dn+ Be) 8-+dy-+49) 
rd+p+h) 


where }+4u+4s > 0, }+4v+4s > 0. Hence, by Mellin’s inversion 
formula (6) and the theory of analytic continuation, 


ain tiy—te-t', Ht») (32") 





ptio 





1 susp (E+3e+ 38) 0 (E+ b+ 38) 
pest Qis+iptiv—-t ads, (3.1 
~a 4 ra+ P49) we 
p—iao 
where p > p+v—A+ 3. This integral is of the same form as (2.3) with 
Diu+iv—t 





x0) ~ TET P+ G+ PW)" 

and 0 < p+v—A+ 3 < 1. Hence, the function on the left-hand side 
of (3.1) may be taken as the function f(x) in the above theorem. 

In 1932 I gave three examples of functions of the type of P(x) in 
the above theorem (9): 

(a) ate DKy,_ (x), 
where X,,,(2) denotes Bessel’s function of the second kind of imaginary 
argument and order m; 

(b) ahi (2); 

(c) IYI (2). 

I now proceed to derive self-reciprocal functions. Taking the 
function (b) given above as our P(x), we have 

g(x) = J (ay)HO-H DT os y(ry)yhet?De-W asty—prau-l BY?) dy 

0 


= gie-n+)) | YY Sys y(2Y) Wy say sav—pru—(3Y") dy. 
0 


Evaluating this integral by a formula given by Varma (12), we have 
tietielt: iw 

1+$u+4hy,1+$4 hr’ 

== 2h-te-tgr+t Fh(1+y, 1+40+)); — 32?) (v > —1). 


g(x) = 20-te-ty 4 2A 
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Using Kummer’s formula 
(a, p32) = &(p—a; p; —2), 
we find that the function 
x’ tte-ie" FL (4A— 40; 1+ 40+ }); 42?) (v > —1) 
is R,. In other words, the function 
gAtin+te-iz" F(—n; 1+-n-+A; $a?) (A+2n > —1), 
is R). 
Deductions from this result: 
(x) Since, when n is a positive integer, 
P(l+a+n) 7p 
(a)(y) — aiaiae v= 
Li? (x) = n'T(a+-1) 4(—n; «+1; 2), 
where the expression on the left-hand side denotes the generalized 
Laguerre polynomial of order n, it follows that the function 
gAt+2ntte—ia* J (n+A( 1x2) 
is R, for positive integral values of n. This result was given by 
Howell (8). 


(8) Since, when n is a positive integer, 


n = xd a 7 a . 
Th(t) = niP+p) i” 1+p;2), 
where 7'ii(x) is Sonine’s polynomial of order n, it follows that the 


function gh+2nthe-da7m (1x2) 


is R, for positive integral values of n. 
(y) If we put n = —A—4, our function becomes 
w-tee" (A+ 45 4s i). 
As a particular case, when A is an integer, this becomes a constant 
multiple of g-te-t"D), (2). (3.2) 
(5) When n = $—A, the function is 
wee", (A—4; $5 42). 
If A is an integer, this is a constant multiple of 
ate" D 4) ,4(2). (3.3) 
The functions (3.2) and (3.3) were given by Mitra (3). 
If we take the function (c) given above as P(x) we get the same 


R) function. 
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4. Example II 
I start with che following form of Mellin-Barnes’s integral (13) for 


Wye m(2): 
Ww. m2) = 


el f P'(s—k)l(—s—m-+})I(— —8+M+4) 6 ae 
P(—k—m+4)l'(—k+m+}) 


Qn 
—io 
where |argz| < 37, and neither of the numbers k-+-m-+} is a positive 
integer or zero. 
From this formula we get 


giutv—lelayyy 
ziuty—Nel W_4-4r-4p-1v,10-p( 42") 
H1—p—) +iw 


2D (E+ 3u+438) (+ vt 38) 0 (3+ 3A—4}s) 2-* ds, 


i BOLE A+ $y) (A+) | 
M1—p-y) iw (4.1) 


where |argz| < #7, and neither of the numbers A+yp, A+v is a 
negative even integer or zero, and |u—v| is not an integer. When 
A = $u—4—1, a pole appears on the imaginary axis. This may be 
avoided by taking the integral along the line (p—ioo, p+-ioo), where 
p>. 
The integral in (4.1) is of the same form as in (2.3) with 
Qi(u+v—s) 
AO) 
PU+P+ iD P+) 

and —l<ptry<l. 


So the function on the left-hand side of (4.1) may be taken as the 
function f(x) in the above theorem. Hence, taking the function (6) 
given in §3 as P(x) and the function on the left-hand side of (4.1) as 


f(x), we get 
g(x) = (ay) HDS ay (zy)yter— ety’ x 


x W4-4)-1p-1y,10-p(dy*) dy 


= ahe-n+)) | PT ry (ry ee W_y_4)-4.-1,40-p(3Y*) Iy- 
0 


Evaluating this integral by a formula given by Erdélyi (4), we get 
glee) = awe 4 ay any n(b2%) 


2 7 
= rity Del Ws ay sray-n(42"), 
3695.10 ad 
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where p+v > —2, v > —1, 4v—}u—A—1 < —}. It follows that 
the function aht2m—tel WT ys ($a), (4.2) 
where A+4m > —1, p > —1, 4m >A, is R). 


This function was given by Bailey (1).* 

If we take the function (c) as our P(x), we get the same R), 
function. 

5. In 1932 I proved (10) that: 

If L(x) is R,, and K(x) is R,, the function 


os x 
a 


ye. - l nf 
Q(x) = | K w)1(") dy [. K(-\Lw) dy 
17 y Io 
0 
is of the form (2.1). 
If we put this function for P(x) in the above theorem, we arrive 
at the result: 


If L(x) is R,, and K(x) is R,, and f(x) is given by (2.3), then the 


function - ; 
% -(xy 
g(x) = [ sw) dy | L(u)K (=) du 
u u 
0 0 
= [ sw) dy | = (4) K(w) a 
7 Ju \u 
is Ry. 0 0 


If in Theorem I we take A = », (2.3) shows that f(x) becomes R,,. 
For example, if in (4.1) we put A = v, the function on the left-hand 
side becomes 


gilh ty—Nele* |v 


ity-ty,4-—p(4 _ sna 42") 
( 


2) — 2h +2m—beiz* YW 
which is the same function as (4.2) with » written for A. 
If we put A = v in (3.1), the left-hand side becomes 


aku ote", +hu—ty, dpe (3%). 


Using the well-known formula 
4 +4,—12 
Wr+t,m(2) = gerie~™, 


we find that the above function becomes the familiar R,, function 
ae bo- ba* 
* Recently, Dhar (3) has given another proof of the fact that the function 
(4.2) is R,. The other Ry function which he gives, 
gh—2m—beket yy j A+3m,m( 32"), 


is the same as (4.2) with ‘—m’ written for ‘m’. 
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Thus, when A = v, Theorem I assumes the simpler form: 
If f(x) is R,,, and P(x) is given by (2.1), the function 
g(x) = [ Play) fly) dy 


2 0 
is R,. 


This theorem was proved by me in 1932 (9). 

It is interesting to compare Theorem I with a similar theorem 
proved recently by the present writer and R. V. Shastry (11). That 
theorem may be thrown into the following form. 


TueorEM II. Jf ° 
P(}+3u-+ $8)P(3+ dv—}s)w(s)a* ds, (5.1) 


where w(s) satisfies (1.2) in the strip (1.3), and 
k+io 
QT (}-+dut+-ds)P ($+ by—4s)P(E-+ 4+ 48) x(8)~* des, 
k—io 
where x(s) satisfies (1.2) in the same strip, then the function 


ao 


g(x) = | P(y) f(xy) dy 


0 


is R). 
If A = v, f(z) becomes R,,, and the theorem assumes the simpler 
form: 
If f(x) is R, and P(x) is given by (5.1), the function 
g(x) = [ Ply)flry) dy 


0 


is R,. 
This theorem was also given by me in 1932 (9). 


6. THeoremMs III-IV. Theorems I-II remain true if w(s) and x(s) 
satisfy the equation 


(8) = —p(1—s) (6.1) 
in place of (1.2). ; 

TueoreMs V-VI. Jf one of the functions w(s), x(s) satisfies (1.2) 
and the other satisfies (6.1) the functions g(x) in Theorems I-I1 are 
_ —R,. 








ON SELF-RECIPROCAL FUNCTIONS 


REFERENCES 

. W.N. Bailey, ‘Some classes of functions which are their own reciprocals 
in the Fourier-Bessel integral transform’: J. of London Math. Soc. 5 (1930). 
258-65, § 6.3. 

——., ‘Self-reciprocal functions involving confluent hypergeometric func- 
tions’: ibid. 13 (1938), 111-12. 

. S. C. Dhar, ‘On certain functions which are self-reciprocal in the Hankel 
transform’: ibid. 14 (1939), 30-2. 

. A. Erdélyi, ‘The Hankel transform of a product of Whittaker’s func- 
tions’: ibid. 13 (1938), 146—54 (2.06). 

. 8S. Goldstein, ‘Operational representations of Whittaker’s confluent 
hypergeometric function and Weber’s parabolic cylinder function’: Proc. 
London Math. Soc. (2) 34 (1932), 103-25 (91). 

. G. H. Hardy, ‘On Mellin’s inversion formula’: Messenger of Math. 47 
(1918), 178-84. 

—— and E. C. Titchmarsh, ‘Self-reciprocal functions’: Quart. J. of Math. 
(Oxford), 1 (1930), 196-231, Theorem 8. 

. W. T. Howell, ‘A note on Laguerre polynomials’: Phil. Mag. (7) 23 
(1937), 807-11. 

. B. Mohan (formerly B. M. Mehrotra), ‘Some theorems on self-reciprocal 
functions’: Proc. London Math. Soc. (2) 34 (1932), 231—40, §§ 4, 8. 

. —— ‘Definite integrals involving self-reciprocal functions’: Bull. Cal 
cutta Math. Soc. 24 (1932), 163-76, Theorem II. 

—— and R. V. Shastry, ‘Theorems connecting different classes of self- 
reciprocal functions’: Proc. National Acad. of Sciences (Allahabad), 7 
(1937), 1-5. 


. R.S. Varma, ‘Some intinite integrals involving Weber’s parabolic cylinder 
functions’: J. Indian Math. Soc., New Series, 3 (1938), 25-33, §§ 2.1, 3.6. 

. E. T. Whittaker and G. N. Watson, Modern Analysis (Cambridge: 4th 
edition), 345, ex. 1. 





ON THE EXPANSION OF A FUNCTION IN A 
GENERALIZED NEUMANN SERIES 


By D. P. BANERJEE (Mymensingh, Bengal) 
[Received 20 July 1938] 


IN expanding an arbitrary function in a generalized Neumann 
series* of products or squares of Bessel functions, it is difficult to 
obtain the expressions for the coefficients in the expansion. In this 
paper I show that the problem can be solved in certain cases with 
the help of the operational calculus. My sincere thanks are due to 
the referee for his kind and helpful criticisms. 
I shall need the following results:f 
f xe-Pe Iya) Iya) de = 50 -# 0d (3) 


(R(p) > 0; R(n) > —1), 


- . nn 
gntle—pz J (Ax) dx = - 


Bp? (Rp) > 0; Rn) > —1), 


0 (2) 
ca g PE | 
ar = Qn > orenet 1 Tale). (3) 


s=0 


I shall also need the following lemma:{ 


Lemma. Let f(z) = > b,,2" be an integral function of exponential 
n=0 


type, so that f(z) = O(e°*!) for some c, and let ¢(w) = > n!b,w-". 


n=0 
Then the series for ¢(w) is uniformly convergent for sufficiently large w, 
i.e. for |w| > M+6 > ‘M >, and, if T is a closed contour surround- 
ing the origin and lying outside the circle |\z| = M, then 


(2) = [ d(w)e™ dw. 


Also, if R(w) > e¢, d(w) = [ f(a)e dz. 
0 
* Watson, Theory of Bessel Functions (Cambridge, 1922), 524. 
+ Gray, Mathews, and MacRobert, Theory of Bessel Functions (Macmillan, 
1922), 33 and 69. 
{ For the proof see Titchmarsh, Theory of Fourier Integrals (Oxford, 1937), 
§ 1.30. 
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Case I. Let f(z) be an even integral function such that |f(z)| << Kee*" 
and suppose that it is to be expanded in a series of the form 
3 a,, J, (Az) I, (uz) 


n=0 


ao 


Let d(w) = [ xev?' f(x) dx (R(w) > ce). 
0 
Then the analytic continuation of this function will be regular for 


w| > M (say) and we shall have, by the lemma, 


J(a*) = a 2e""d(w) dw 


wr-lewe-W— pw Qype — Awd (w) dw 


r 


] ae [A 
| w-lewe—Or- p/w FP) dv, 
271 J 2w 
r 


3 A 80)),) 
where Piz) = Eel H weg), 
, shee ote Fle) j . ee Au| 
But wd(w) is analytic for |w| > M; thus F(z) is analytic for |z| <=_, 


and within this region . 
Al 
F(z) _— > Ay J,,(z), 
n=0 

where* 

<(in] 7, ‘ 

(n—2m) apn 

I (9) 5n- om (” 


ae ! 
i=, (n—2m): 


ay, = F (0), =e 9 
the series being uniformly convergent for 


r(3*) S a, (5 1, 
2w 2w 


n a, 


It follows that 


the series converging uniformly on [. Since w-te*-™’-r)4w js 
bounded on I’, we van substitute for F(Aw/2w) in (4) and integrate 
term by term, giving 


x 
f(x) — > Sy a evr (A? a Au dw 
2m 2w] w 
0 “ 
I 


n= 


=> a, T(ux'), 


* Watson, loc. cit., § aa (3). + Gray, loc. cit., p. 35. 
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the last result being obtained from equation (1) by applying the 
lemma to the integral function J, (Az?)I, (ux!). Hence we have finally 


=> a, J, (Ax)I, px), 


where a, is given by (7). 
. Let f(x) = x"J,{(A2—p?)!2}. Then 


Example 1 


‘te [ e-watynt1 {(\2—p?)bx} dar 
0 
so (A2—x2)"/2( 2a) —"—“le-O*— #4, 


by (2). Substituting into (5) we have 
2 ,2\ nl2 
r= Cae) * 
p2\nl? < (n+ 2s) at 
sai (ae 2 )" bs en Sovalth 


=0 
by (3), and thus it follows from (6) and (7) that 


2" J,f(X?—p2) ha} 
-_ n[2 © 9 
a -) pc eine INF seal Ty soe(uex). (8) 





A?y? s! 


Example 2. Putting ip pi p in (8), we have 


or e+ p2)h2} 
? nl2 = —1)*(n+2s)(n+s—1 
=e re seas). 

a8 — s! 
Case II. Suppose that f(x) satisfies the same conditions as before 
and that it is to be expanded in the form > a,,2"J,,(x). We have 





n=0 


] 
f(x) = — | 2e""d(w) dw 
2Q7i 

i 
Lf 4 -teme-ttedangtite 
. | w-lewr-Ue2yelMud(w) dw 
272 | 
“ 
Lf gteus-titw 19 
— | w-lewr-liw F(1/2w) dw, 
2m | 
rs 


where F(z) = z-1e!* (3). 
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Since F(z) is analytic for sufficiently small z, we have 


i a) 
pe — FM (0 
f(x) = S —— | etr—4w(2w)-"-1 dw 
La 2 ml ott 

n=0 rr 


. 


— = : 
_ << ered (a! 
> am, (al), 


n=0 


by the formula reciprocal to (2) (with A == 1). Hence 


n! . 


. = nn) 
j(x) = 2 x" J (x). 
=() 
The analysis is easily justified as in case I. 
Example 3. Let f(x) = 2". Then d(w) = 3m! w-"-! and 
J\ 


F(z) = m!(2z)me? = 


x 


It follows that 2m = Om > ( — (5) J},(2). 
n—m)!\2 


n=m 
Case III. When case II fails, or when the evaluation of 4(w) is 
difficult, we may sometimes proceed as follows. Let f(z) be an even 
integral function such that |f(z)| < Ke°*', and let ¢(w) be defined 
as in the lemma. Since b, = 0 (n odd), wd(w) is a function of w?. 
We have » 


‘ | e”“d(w) dw 


. 


oi | (l = (1+-w?)!d(w) dw 
f 


——, ! = —— fF : -| dw, 
J (1l+w?)! \l+w? 
s 


where ; r *)'b(w) ( a i wd(w). 


rhis is an analytic function of ; ; for sufficiently large w; hence 


oe ‘(n ¥ wr 
fle) ‘= i 0) a | , ¢ - dw 


nm! 2n0 J (1-+w?)"4 
n=0 V 


ive) 


> (7) (()) ar 


av nm! 2"P(n+4) 
n=() 7 


x" J, (x), 
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2"T'(n+4) 
(1+ w2yn 
Example 4. If f(x) = cosx, then ¢(w) = page and 


since &”"J, (x) and are Laplace transforms.* 


1 . .(2n—3) 
Fo) — aa 2a at (1+-w?)” “ 


— ()"d, 
H OF == Bas a 
ence cos x 2a '(2n_-1)° 
Example 5. If f(x) = xsina, then d(w) = 2w/(1+-w?)? and 
1 2Qw = 1, (2n—3) 
Wal = (i+w) an 2st 2 2" iol 


\ n\(4n?—1) 


* Titchmarsh, loc. cit., p. 207. 


Hence reins = fad (x)—2 S er n+(% =} 








THREE TRIPLE INTEGRALS 
By G. N. WATSON (Birmingham) 
[Received 3 June 1939] 


The desirability of investigating the triple integrals 


I dudvdw 
1 = 3 
1 — COS U COS V COS W 
000 
TW 7 
l [| [ dudvdw 
- 7 ? 
7 | J J 3—cosv cos w—COSs Ww COS U— COS U COS UV 
000 
TW WT 
re l | | j dudvdw 
5 3 J J J 3—cos u—cos v—cos w 
000 


has arisen as a consequence of their having appeared in a recent 
paper* in ferromagnetic anisotropy by F. van Peype, a pupil of 
H. A. Kramers. The problem of evaluating them was proposed by 
Kramers to R. H. Fowler who communicated it to G. H. Hardy. The 


problem then became common knowledge first in Cambridge and 
subsequently in Oxford, whence it made the journey to Birmingham 


without difficulty. 

It is possible to express all three of the integrals in terms of simple 
surds, the number z, and certain complete elliptic integrals. The 
elliptic integrals which occur in J, and J, are easily expressible in 
terms of gamma functions whose arguments are simple fractions. 
This is not the case, so far as I know, with the elliptic integral 
occurring in J,; though this elliptic integral is quite familiar, since 
its modulus is one of the simplest of the set of singular moduli. 

The value of J, appears to be fairly well known, but, for the sake 
of completeness, I give two methods of determining it. So far as I 
know, no method of evaluating J, and J, has hitherto been discovered ; 
my methods involve, in the one case, the use of a certain expansion 
for the square of a complete elliptic integral and, in the other case, 
the use of a formula discovered by Bailey connecting the product 


* Physica, 5 (1938), 465. 
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of two hypergeometric functions with a hypergeometric function of 
two variables of Appell’s fourth type. 
The results which I have obtained are as follows: 


_4K3_ Pq) 


> 
7 473 


ql 


os K2v3 3184) 
_ . 

I, = (18+12V2—10v3—7+6)(2K,/z)?, 
where K,, K,, K, denote the complete elliptic integrals with moduli 


sin 45°, sin 15°, (2—~v3)(v3—~2) 


respectively. From these results it is easy to compute the numerical 
values of J,, J, and J;, namely 


I, = 1-39320 39297, 
I, = 0-44822 03944, 
I, = 0-50546 20197. 

I start with all three integrals in the same way, by writing 


tan ju = 2, tan dv = y, tan jw = z, 


and then, regarding x, y, z as Cartesian coordinates, I change to polar 
coordinates by the transformation 


x = rsin@ cos ¢, y = rsin@sin 4, z= rcos6; 


after writing 26 = #, I next bisect the ~-range of integration; the 
integrands at corresponding points in the two parts of the range 
being equal, this procedure gives 


in : iv 
| (...)dé = t| (...)dyb = | (...) dep. 
0 0 0 

There is no difficulty in justifying the various changes which are 
made in the order of carrying out the integrations, because all the 
integrands are positive throughout and the triple integrals are trans- 
formed into integrals which are evidently convergent. 
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By treating J, in the manner just described, we get 


dadydz 


pes 
a” 7 | | | (1+-a)(1+-y?)(1-+2?) — (1 —a®)(1—y?)(1 —2?) 


000 


x x x 
4ff dadydz 
~ % 3 2 S| 2,202 
mw Jy) e+ y®+2?+ aty22* 


00 0 


hor har x 
sin 0 drdédd 
1+r4 sin‘ cos?6 sin?¢ cos" 


4 “- a7 a7 
a7 | | | 
0 0 0 


oS e sin 6 drdédis 


a J J J 1+4rsin‘6 cos?6 sin®s 
0 0 0 


In the innermost integral replace the variable r by a new variable 


¢t defined by the formula 
t = rsin 6,/(4 cos @ sin #); 


this transformation gives 
dtdédis 


ae... ; 
1 ‘ n /($ cos 6 sin yb) 


a | dt [ aor dis 


7 1+¢ J Vvcosé J vsing 
0 0 0 
Now it is an easy exercise in contour integration to prove that 


x 


dt 7 
1+7#4 2/2’ 


0 


and it follows from known properties of the first Eulerian integral 


that ho hor 
. we [ dis 1yr(4 P?(4) 

V cos 6 V sin 2T(3 2,/(27) 

0 0 
By combining these results we immediately get 
rc 
Pe ee, 
4773 


which gives one of the values stated for J,. 
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The second method of dealing with J, is to effect one of the inte- 
grations forthwith and get 
TT 
he dvdw 
n° JQ (1—cos*v cos*w)’ 
We now write 
1.3.5...(2n—1) 
a, = 1, a, = — —— 
2.4.6...(2n) 
for brevity, expand the integrand in ascending powers of cos*v cos*w 
and integrate term by term; we thus get 


(n > 1) 


1 S 
i= wD, er cos*"w dvdw = > ae 


the term-by-term sisal being justified by the convergence of 


the series on the right. 
Now it is known* that, when 0 < k < sin 45°, 


K? = Yn? 3 ad(2kk')"; 
n=0 


and so, taking k = sin 45°, we find that 
4K5 


=F S 


i= 


which gives the other value stated for J. 

The numerical value of J, is now immediately derivable from 
results contained in Legendre’s tables of elliptic integrals and 
gamma-functions. 

We now turn to J,. When we transform it in the prescribed 
manner, we get 


«oo 


sf dadydz 
=| 
00 


3 TT (1 +2%)— ¥ (1—y*)(1—2*)(1 +24) 





iL= 
9 ow e®w —— dxdydz 
m1 7 ay? 22+ y222+ 22224 ay? 
a a) 


2 ¢ sin 0 drdddd 
7 | | | 1+ -r? sin?0(cos?6+ sin?6 sin*¢ cos*¢) 
0 


0 0 





* Expansions of this kind have been discussed by various writers ; I dealt 
with them myself many years ago at some length, Quart. J. of Math. 39 (1908), 
27-51. 
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sin 6 drdédys 
J J 1+7?sin?0(cos?6-+ } sin?@ sin) 
dds 


cos*6-+- } sin? sin*/) 


dtd 
mJ J 1 +8)1+ fe sins) 
where, at the last step, tan@ has been replaced by t. 

We have now to introduce complete elliptic integrals of the first 
kind in which it is convenient to make explicit mention of the 
modulus. We therefore write K(k‘) to denote the complete elliptic 
integral of the first kind with modulus k and we write K’(k) to denote 
the complete elliptic integral of the first kind in which the modulus 
k’ is the modulus complementary to k, so that K’(k) = K(k’). 

se(u,k) 
dt 


Since - <r. 
MO+8)(1+k)} 


x 


At 
‘e he TA [ a K k — vr k’ : 
we have LL E\L-Ek?P)) X(h) («) 


© 


0 
whence we obtain the preliminary result 


}7r 


L=— [ K’(4 sind) dy. 
a 


ry 


0 


Now, with the meaning already assigned to a,, we have 
f=) * e nt 
=8 
K'(k) = ¥ adfp(n+1)—p(n-+})—log kk"; 
n=0 


and, for values of k between 0 and 1, every term in the series on the 
right is positive. We consequently have 


‘te 3 ; 

. A az | {y(n + 1)—d(n+ $)—log(4 sin %)}(4 sin fs)?” dy, 
mr n=0 b ‘ 
0 


term-by-term integration being permissible because the series on the 
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right proves to be convergent. To evaluate the general term on 
the right, we observe that 

2 T(in+$+e) 

a2s = 1 = kik2n wide: ees inte 
aniy(n+ ) (n+) Og ks lap T(n+1+e) )I.. 
and so we have 


in 


a2 | {yb(n+1)—o(n+ 4) —log($ sin #)}(} sin ys)?" dip 


__ ay ole (nests f (4sin #) )an+e a) 


Nir (n+1+ & 


de \22"+€+1P (n+ 1+e)P(n+1+ fe) 


= }7a3{3yh(n+ 1)—3yb(n+ 4) +2 log 2}($), 
the interchange of the order of differentiation and integration pre- 


senting no particular theoretical difficulty. We have consequently 
arrived at the result 


NTT 


31 2(5 T(n+ +1 ileal MI... 


hot > atten) Hon) + 20g 21 


Now in my previously cited paper it is proved that, for 
0<k< sin45° 


(and by symmetry also for sin 45° < k <1), we have 
K(R)K(R’) = dS ah{ (n+ 1) —h(n-+ 3) —2 log (2KK")}(2KK)™ 


When we compare this expansion with the formula just obtained 
for J,, it is clear that 
K (ky) K (ky) 


2 


L= 

7 

where k, is the modulus between 0 and sin 45° which satisfies the 

condition 2k,k; = 4, so that k, = sin15° and k; = sin75°. If we 

now use Legendre’s well-known result that A(sin 75°) = K(sin 15°)v3 

together with his formula expressing K(sin 15°) in terms of gamma- 
functions, we immediately find that 

3S K%(sin15°)v3_ —- 3P8(3). 


= =. l ? 
red 914/34 


these are the formulae for J, stated at the beginning of this paper, 
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and the computation of the numerical value of J, is merely a matter 
of employing a table of gamma-functions or of complete elliptic 
integrals. 

We finally consider J. When we transform it in the prescribed 
manner, we get 


aoa 
a 


ah {J J dadydz 
a JJ > 2742 > y22?+ 32%y22? 
0 


sin 6 drdéd¢d 
1+ 2r? sin?6(cos*@+-sin?6 sin*¢ cos*¢) +- 
+ 3r4 sin*6 cos?6 sin*¢ cos*d 


r sin 6 drdédys 
| | | 1 + 2r? sin?0(cos?6+- } sin20 sin’) + 3r4 sin4@ cos?0 sin®y" 
0 0 0 
Now replace r by a new variable ¢ defined by the equation 
rsin@ = tv2 and then perform the integration with respect to @. 
This procedure gives 
br in 
eo “1, iy [ { dtdods 
. 1 +#?(4 cos?6-+- sin?6 sin*ys) + 3¢4 cos?@ sin*s 
0 0 0 


7 x 
+ | dtd 
q {3 (1+ 442+ 3¢4 sin*ys)(1-+-é? sin®y)} ° 


0 0 
Next take tan, as a new variable é, and then write € = ¢/,/(1+1#°) 
and perform the integration with respect to t. We thus get 


dédt 


I; f(1+ €#-+. 442+ 47? i 3t*é?)(1 + €*+ t7€?)} 


dtdt 
14424 24 32f)(1+2)} 


x(t 
J 44-32?) /{(1+-0)(4+30")} 


- 1 dy 
- |; . j . = 49 
V(4—sin?x)/ ,/(4—sin*x) 
0 


by writing ¢ = tan x. 
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We now replace K’(k) by a differential coefficient, as in the corre- 
sponding discussion of J,. The expression which it is now convenient 
to take is not the one used in dealing with J, but 


eat 
-) — te: oth 


The changes of order of integration, differentiation, and summation 
are even more easy to justify than in the case of J,, since the range 
of values of k now under consideration extends only from } to 
1/V¥3 and it has not now got 0 or 1 for an end-point. We thus 
find that 


ia 
bk Minti 1 
ait -3iz (2 l'(n+1-+ 2e) aaa) | 
= bs 


vafdS/rmttte fay 
i. Re: n! TP(n+1+: 20 J w 4—sin®y)"+i+« ey 


The next step is the integration of the various powers of 4—sin*y. 
We define c to be the smaller root of the quadratic equation 


c2+1= 4c, 
so that ¢c = 7—4v3, and we then have 
4c(4—sin*y) = 2¢(7+-cos 2y) 
= 1+2ccos2y+c? 
= (1+-ce%x)(1+4-ce-2x), 


whence it follows that 
1 ho (4c)" ttt 
(4—sin?y)"+#+© — {(1+-ce?*x)(1+-ce-2x) In tite® 


We now expand each factor in the denominator on the right in 
a series of ascending powers of c, multiply the two series together, 
and integrate. Since the integrals of the products of unequal powers 
of e?*X and e-*** are all zero, this procedure gives immediately 


47 


dx 
(4—sin?y)"+!+€ 
) 


= }n(4c)"+#+* FF (n+3+e,n+4$+¢; 1;c?%), 


3695.10 T 
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whence it follows that 


> (Bet 


n!T(n+1 


sos 
n=0 


+ 2e, 1; 4c, c?), 


where §, denotes the fourth kind of Appell’s hypergeometric func- 


tion of two variables. 
We now avail ourselves of Bailey’s result* that 
© . Q . 
4 x, Bs y, x-+-P y+1;a(1—b), b(1 a)) 
= oF (a, By; a)oFy(a, B; a+B—y+1;), 
which, for a given value of the (complex) number a whose modulus 
is less than unity, is valid for values of b in the domain surrounding 
the origin for which the condition 
‘y |a(1—b)|+./ |b(1—a) 
is satisfied. When 
a(l- 
> 
a, 
we have  |a(1—b)|+ 6v3 < 0, 
and so the requisite condition is satisfied. Moreover, when a has 
any given positive value between 0 and 1, the domain of values of 
b surrounding the origin within which the condition is satisfied con- 
tains that portion of the real axis which extends from the origin to 
1—a. We therefore must select as the values of a and 6 that pair 
of solutions of the equations 
a(1—b) = 28—16v3, 
b(1—a) = 97—56v3 
which satisfy the condition a+b < 1, assuming such a pair to exist. 
There is, however, no difficulty in verifying that 
a = 24v2—14v6+ 20vV3— 34, 
b = 24V2—14V6— 20vV3+-35 
are solutions of the simultaneous equations such that 


a+b—1 = v4608— v4704 < 0. 


* W.N. Bailey, Quart. J. of Math. (Oxford), 4 (1933), 305-8. 








THREE TRIPLE INTEGRALS 
We now write 
k= = (2—~3)*(v¥3—v2)?, 
k2 = 1—a = (2—v3)9(v3+ v2), 
and we immediately find that 
. é (Maree 


id Lie Lite: 14%: ki? 
de\ Tata) Pate bes 14265 bs) x 


m2 
x F($+e, $+; 1:19) : 
e=0 
An application of Euler’s formula to the second hypergeometric 
function on the right gives 
F(4+e, $+; 1; 3) = (1—K)-?-*F(4—e, +; 1; —K3/k,?); 
now F(4}—e,4+e«;1; —k3/k,?) is an even function of «, and so its 
differential coefficient with respect to « obviously vanishes with e. 
Consequently, when we differentiate the product of the functions 
F(4—e, $+¢; 1; —k3/k,?), 
T'?(4-+-e)kg!+¢ o 
=. : F(4-+e, 4+; 1+ 2e; k,?), 
I'(1+ 2e) Gre: 3) 
we immediately obtain the formula 


1; —k3/k,?) x 


7 d [(4+-e)k3! +2€ 
de\ T(1+2e) 
A second application of Euler’s formula now gives 


I, = _h% | F(3, 35 1; k}) x 


rv 2 
2(4 €)ky?* Wi} 1 ¢ , 
— ——— ste, 5+ ;1+ 2 > ke? 
x [a(- P'(1+2e) Fite ite Ks’) 


_ 2ky ky Ky Ky V2 


2 > 


F(4+e,4+¢; 1+263k3))| 


€=0 


e=0 


7 


where K, and Ky, naturally denote the complete elliptic integrals with 
moduli k, and ks. 

Now, if 7,, denotes the parameter associated with the modulus f,,, 
that is to say that 7,, is the quotient 
= ik nl Ky, 


Tn 
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where Kj, denotes the complete elliptic integral whose modulus k;, 
is complementary to k,,, it is known* that k, and k, are the moduli 


for which 


* 19 


TT? = iv6, Tz = W/§- 
The quarter periods K, and K, are consequently connected by 
Landen’s transformation, so that we have 
Ky = (1+k,) Ko, 
in addition to the previous relation 
K; = K,,/3. 
By means of these results the formula for J, reduces to 


2ke ke + vo 3 e 
hw 2h, ko(1 ft V3. 


> > 


ae 
and, when we replace the moduli by their values in surd form, we 
immediately find that 
I, = (18+12V2—10V3—7+6)(2K,/7)?, 

which is the value stated at the beginning of the paper. In com- 
puting the numerical value of J; from this result, I found it quickest 
to use the formula 

2K,/ma = (1+ 2e-7%8 + Qe-47"6 4 Je-Smv61 2 
for the preliminary computation of 2K,/7 = 1-00182 06717 32. 


* See e.g. Greenhill, Proc. London Math. Soc. (1), 19 (1889), 351. 








ON TCHEBYCHEFF POLYNOMIALS 
By D. R. DICKINSON (Cambridge) 
[Received 10 June 1939] 


In §§ 1-4 I shall reproduce some definitions and some known 
results concerning 7'-polynomials. 


1. We say that a system 
po(), ,(2),..., $n(2) 


of n+1 real functions, continuous in the closed interval (a, b), form a 
7T-system of order n in (a,b) if no linear combination 


P(x Ao $o(2) «)+A, (x )+.. +A, $,(2), 
with real a has more than v distinct zeros in this interval 
unless the coefficients are all zero, and we then call P,(x) a T'-poly- 
nomial (or simply a polynomial) of the system. 

Let @ be a zero of P,(x). Then we say that is a simple zero of 
P. (x), if P,(x) changes-sign as x passes through %. If P,(x) preserves 
the same sign as x passes through %, we say that # is a double zero 
of P,(x). If P,(#) has a zero at an end-point of the interval (a, b), this 
zero is, of course, counted as simple. In this paper I am concerned 
with the problem of constructing 7'-polynomials with previously 
assigned zeros. 

2. Let 2%, X},..., x, be any n+ 1 points of (a,b). We write 


D, (Xo, %;-- pret Po(%) Gil%) - - + Pn(%)|- 
$o(2) ita - + © Gul®%) 


dole) dy Ly)» + + PnlFn) 
THEOREM 1. A necessary and sufficient condition for $o(x), $,(%),..., 
¢,,(x) to form a T'-system is that 
DA dog, By 5.15 By) F O 
provided only that x, x,,..., X, are all distinct. (1) 
Suppose there exist n+ 1 distinct points such that 
D, (Sq, ¥q+++y%_) = 0. 
Then, by a theorem on determinants, there exist numbers Ag, 
A,,, not all zero, such that 


s A; $,(x) = 0 
Ko 








278 D. R. DICKINSON 
for x = 2p, X,..., %,- But this is impossible for a 7'-system. Hence 
the condition (1) is a necessary condition. 

On the other hand, suppose that condition (1) is satisfied. Then, 
ET Bq, Za500es x, are all distinct, the only solution of the equations 


is A, = © & = 0, ]...., n). Hence our functions form a 7'’-system, 
and the condition is a sufficient one. 

3. We shall suppose from now on that the functions {¢,(7)} form 
a 7'-system in the interval (a,b). Now let xo, 2,,..., x, be n+1 distinct 
points of (a,b). Then D, (x, x,,...,%,) is of constant sign provided that 

6% % < ..< 2, &8 

and changes sign if we interchange the inequality between two 
consecutive x;. It follows that the 7'-polynomial D, (x, 2, %9,...,%») 
has n simple zeros at the points 2 = 2, %y,..., x,. Further, any 
T'-polynomial having zeros at these » points is a constant multiple of 
DAS, B55 Bigg s+-3 Bn) 

4. THEOREM 2. If 8 is the number of simple zeros and d the number 


of double zeros of a T'-polynomial P, (x) of our system, then 


st+2d < n. 
Let 2), Xo,..., % be the simple zeros of P,(x) and y,, Yo,..., Yq its 
double zeros. The sign of P,(x) in the neighbourhood of each y; is 
constant and so, taking n—s distinct points x}, x5,..., #7, lying to the 


left of all the y;, we can construct a 7'-polynomial 


Q,, (x) SL Ck 4 ae ae oe ee ® 


— “ - 


having » simple zeros at the points 


, 


ye eee ee ee 
whose sign in the neighbourhood of each point y; coincides with that 
of P,(x). Denote by 26 the least distance between two zeros of P,(x) 
and write 
M; bound {|P,(%)|} (yj; < # < y;+8), 


M; = bound {|P,(x)|} (y;j—-8 < # < y;), 
M = min(M;, M;) (j = 1, 3...., d), 


N = bound {|Q,(x)|} (@a<a <b). 
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Choose a positive number A such that AV < M. Then to any point 
y; there corresponds at least one zero of the polynomial P,,(x)—AQ,,(x) 
in each of the open intervals (y;—4,y;) and (y;,y;+6). Apart from 
these zeros the points 2,, 2 ,..., %, are zeros of P,(x)—AQ,,(x). Hence 
the polynomial P, (#)—AQ,,(x) has at least s+ 2d distinct zeros, and 


therefore s+2d <n, 
as was to be proved. 

5. We now come to the main problem of the paper. 

THEOREM 3. Jf s+2d =n, and we are given s+d distinct points 
Ly Loyeeey Les Yr» Yos---> Yq SuCh that 


(¢ = 1, 2...., 8), a<y,<b (j=1, 2....,d), 


then we can construct a T-polynomial of our system having simple zeros 
at the points 2, a,..., x, and double zeros at the points y,, Yo,---, Ya-* 


First suppose d = 1. Take points f,, y, such that 


B<%< "1% 


and, further, such that the interval (f,,y,) contains none other of 
the given points. We can construct a unique 7'-polynomial taking 
the value unity at the point 8, and having n simple zeros at the 
points y1, 24, X,..., V, Y¥,; Denote this polynomial by P(x), and 
denote by Q(x) the unique 7'-polynomial taking the value unity at 
the point y, and having simple zeros at the points f,, 21, %9,..-, Zs Y1- 
Observe that -the polynomials P(x), Q(x) are of different sign in 
the interval (8,,y,) and have the same sign outside this interval. 


Now let B,<y¥<y Or yW<y¥<y 


Qs?(y) 


and write Ay) = — P2y)’ 


Then A,(y) is non-negative in the ranges of y considered. Further, 
the 7'-polynomial Q0)(x)+A,(y)P@(x) takes the value unity at the 
point y, and has n zeros at the points x, %9,..., %,, Y,, y- It follows 
that A,(y) cannot assume the same value for different values of y. 


* This result has been given by 8S. N. Bernstein, Eatremal Properties of 
Polynomials (in Russian) [1937], p. 9, lemma II. Bernstein’s proof, however, 
does not seem to be complete, and I give here a proof on entirely different 


lines. 
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Hence, since A,(y) is continuous, it is monotone increasing (in the 
strict sense) in the ranges of y considered, and the limits 
A, = lim A,(y), A, = lim A,(y) 
ae y>"—0 yr"n+0 

exist finitely and 0<A, <,. 

Choose yp, so that A, < x, < A, and consider the 7'-polynomial 

ROY (x) = Q(x) +1, PH (2). 

R® (a) has a zero at the point y, and, since the function A,(y) never 
assumes the value j,, it can have no other zero in the interval (8,, y,). 
Further, R(x) takes the value unity at the point y, and the positive 
value p, at the point B,. Hence y, is a double zero of R(x). But 
R® (x) has also zeros at the points > Xpy,..., &, and—by Theorem 2, 
or from the fact that P(x), Q(x) are ae the same sign outside 
the interval (8,,y,)—these must be post’ zeros. R(x) therefore 
satisfies the required conditions. 

Observe that the solution to the problem is unique (except for 
a constant multiple) if and only if A, = 4,. We shall return to this 
point later and show by an example that this is not always the case. 

We now pass to the case d = 2. Take points f,, y. such that 

Bs < ¥2 < Ye 

and, further, such that the interval (f,,y,.) contains none other of 
the given points. By the preceding argument we can construct a 
T-polynomial P‘?’(x) taking the value unity at the point B, and 
having a double zero at the point y, and simple zeros at the points 
Yo, ©, Lp,..., Vy, Yo. Similarly we can construct a polynomial Q{?)(x) 
which takes the value unity at the point y, and has a double zero 
at the point y, and simple zeros at the points By, 21, %9,..., Yo. 


— Wy 
ps 2(y 
P(x) and Q!?)(x) are of different sign in the interval (8,, y.) and of 
the same sign outside this interval. It follows that A,(y) is non- 
negative and that the polynomial Q‘(x)+-A,(y)P‘ (x) has a double 
zero at the point y, and simple zeros at the points 2, %9,..., Xs, Yo; Y- 
Hence, by Theorem 2, A,(y) cannot assume the same value for 
different values of y and we deduce, exactly as before, that the limits 
A, = lim A,(y), A, = lim A,(y) 


Yy>¥2—0 YU. t+ 0 


Now let Boa Y¥< Ye OT Yeo<Y¥< Ye 
and write A(y) = a 
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exist finitely and 0<A,<A,, 


and that, if we take any number p, such that A, < p, < Ag, then 
the 7'-polynomial 
RR) = WP@)+mPP@) 

satisfies all our requirements. 

The general case of the theorem follows immediately by repeating 
the above argument. 

6. Finally, we give an example which shows that the solution to 
Theorem 3 is not necessarily unique. 

THEOREM 4. The functions 

1, |x], sin x 

form a T-system of order 2 in the interval (— 42, $x). 

We have to show that, if a, b, c are not all zero, then the function 

f(x) = a+b\x|+csinzx 

cannot have more than two zeros in the given interval. This is 


evidently true if c = 0. Suppose, then, that c is different from zero. 
Without loss of generality we may suppose c = 1. Then 

f(x) = b+ cosx (x > 0), 

f'(—9) - —b+1, 

f(x) = —b+cosx (x < 0). 
Suppose b> 0. Then, if b <1, f(x) is strictly decreasing in the 
interval (— 47, —arccosb) and strictly increasing in (—arccos b, $7); 
while, if b >1, f(x) is strictly decreasing in (—}7,0) and strictly 
increasing in (0, $7). Similar arguments apply if b < 0, and in either 
case f(x) cannot have more than two zeros in the interval (— 4, $7). 
If 6 = 0, then f(x) can have at most one zero in this interval. This 


proves the theorem. 
Now write : 


P(x) = = |x|—Jsinz, 


Q(x) = “ |x|+4sinz. 
We have P(—42) oS 1, P(0) a P(4n) te 0; 
Q(37) = 1, Q(0) = Q(—}7) = 0. 
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- Z 
A = lim _ ey) Ses 
~ poo Ply) 


i=in — 2 _* 
y>+o Py) 


Hence, if we take any number p such that 


and write R(x) = Q(x) +pP(2) 


_ : (1+p)\x|+3(1—p)sin a, 


7 


R(x) will take the value unity at the point $7 and have a double 


zero at the origin. 











FREE STREAM-LINE FLOW PAST A VORTEX 


By N. SIMMONS (London) 


[Received 16 February 1939] 

“1. Introduction 

THE determination of the influence of boundaries on the field of 
motion around a vortex is a problem which assumes an important 
aspect in wind-tunnel theory. This aspect has interested a number 
of writers*, and solutions have been obtained for both two-dimen- 
sional and three-dimensional fields. The methods used are, however, 
open to certain objections. Of these the chief relates to the use of 
image systems. Where the boundaries are fixed walls, on which the 
vanishing of the normal component of velocity is demanded, the pro- 
blem may be resolved by this means. Where, however, the boundaries 
consist of free stream-lines, on which the condition imposed is the 
constancy of the pressure, the image method, strictly speaking, is no 
longer justifiable. Nevertheless, previous writers have used image 
systems for this case (of a type different from that for fixed 
boundaries), taking them to be permissible when the squares of the 
‘added velocities’ may be neglected, and the deviations of the 


boundary from a straight line are small.+ If the boundaries consist 
of both fixed walls and free stream-lines, as in the case of a jet, or 
a wind tunnel with open working section, the justification is less 


obvious. 

In the present paper I solve the problem of two-dimensional flow 
round a vortex, the fluid issuing from between two parallel walls, 
making no assumptions beyond those usual in free stream-line theory. 
The problem of a stream bounded by free stream-lines only, which 
can be deduced as a limiting case from that of the jet, is also treated 
directly. The method used is not entirely new, but is similar to that 
employed in problems of free stream-line flow with singularities by 
Hopkinson (2), which paper seems to have attracted little attention. 

It is proposed in a subsequent communication to deal with im- 
portant particular cases, to show the application to wind-tunnel 
theory, and to compare the results with those of L. Prandtl. Also 
deferred to a later stage is the discussion of the case, approximating 
more closely to a wind tunnel of open working-section, where the 


* See Durand (1), 2, 238 for a detailed list. + Durand (1), 2, 241. 











284 N. SIMMONS 


stream issues from the orifice formed by two parallel walls and 
returns to a collecting channel of similar form. 


2. Stream issuing from between two parallel walls: pre- 
liminaries 
Fluid issues as a two-dimensional jet from between two parallel 
walls CB, C’D (Fig. 1), the free stream-lines being BA, DA’: the 
points C, C’, A, A’ are all at infinity. There is a vortex in the jet, 
as shown. 











A 
a B 1 
hy 

=. Con 
velocity z-plane ae 
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Fic. | 
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For convenience, we take the distance between the walls as z units, 
and the initial velocity of the stream as unit velocity. This is tanta- 
mount to choosing special units of length and time, and so causes 
no real loss of generality. Let the coefficient of contraction of the 
jet be c: then the free stream-line velocity is 1/c. Let the strength 
of the vortex be 2«z. 

Take for origin O the mid-point of BD, with Ox in the direction 
of the undisturbed stream, and Oy perpendicular. Let z denote the 
complex variable in this plane. B, D are the points + (d+ 41), where 
2d is the overlap of the two walls. 

Let Q denote the magnitude of the fluid-velocity at any point, 
and « its inclination to Ox. Let w represent the complex potential, 
on that convention in which 


Qe-' — 


dw 
dz° 


(2.1) 
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We seek a conformal transformation « = u(z) which shall repre- 
sent the area of motion in the z-plane on an area in the u-plane in 
such a way that the boundary ABCC’DA’ may transform into the 
real axis. This transformation can always be chosen so that the fixed 
boundaries correspond to the segment of the real axis joining +1, 
while the point at infinity to which the stream proceeds (A, A’) 
corresponds to u = 00. The free stream-lines transform to the por- 
tions of the real axis of wu outside the segment (—1,1). Let (C,C’), 
the point at infinity whence the stream comes, correspond to Uy 
(Fig. 2). 

The singularity in w at the point z, representing the vortex will 
transform into a similar singularity in w at u,, the point into which 
z, transforms. Motion in the u-plane represented by the complex 
potential w is therefore due to a vortex at u, of strength 2«m7 and 
a source at u, discharging fluid at rate 7, while the real axis is a 
bounding stream-line. Taking account of the image system we obtain 


w = log(u—uy)—tx log(u—u,)+7%« log(u—a,)+-constant, (2.2) 


ai dw _ “eS tx(o,—4)) ba quadratic __. (2.8) 
du U—Uy (u—U,)(U—%,) — (U—Uy)(u—U,)(U—%,) 





The quadratic numerator has real coefficients, and so its zeros are 
either conjugate complex numbers, or both real: we shall say that 



































X denotes stagnation point 
Fic. 3 


the motion is of the first or second type accordingly. A zero corre- 
sponds to a stagnation point if it falls in the upper half-plane. Thus, 
in motions of the first type, there is one stagnation point which 
occurs at an internal point of the fluid (Fig. 3). In motions of the 
second type, there are two stagnation points, which occur on the 
boundary of the area of motion (Fig. 4). These must clearly be 
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situated on the fixed boundaries, and, considering the changes in 
the direction of flow along the walls, we easily find that both stagna- 
tion points must lie on the same wall. This is the upper or the lower 
according as x = 0, the zeros of the quadratic lying, accordingly, in 


the segment (—1, uw») or (wu, 1). Figs. 3 and 4 are drawn for the case 


x > @. 











Fia. 4 
The first type is of the greater interest and will be developed fully. 
The method of dealing with the second type is similar and presents 
nothing novel. The condition for the motion to be of the first type 
is found to reduce to 
(1—K?)S{u,}+ 2«[ R{u,}—uy] > 0. (2.4) 


3. Development of the method 
Define a function Q by the relation 


ie dene. (3.1) 


Then, from (2.1), Q = —log Q+ia. (3.2) 

On the real axis, between +1, 2 must be real, while, outside this 
segment, the real part of Q must be constant. In the neighbourhood 
of wu = +1, 2 must be expansible in the form (+ 1—vw)! x regular 
function. At infinity Q must be finite. If uw; denote the stagnation 
point, Q must have logarithmic singularities at u,, uj, with coeffi- 
cients of the logarithmic parts +1 respectively, so that the relation 
between z and w may be conformal at these points.* Hence the 


function 


i — 
v ‘du 


| 


- / -~/ 
Uu—U, u—t, uU—Uy Uu—Uy, 


dQ poe VU—u’) a) 


is regular at every point of the upper half-plane, and real and con- 
tinuous on the real axis. It can therefore, by the Riemann-Schwarz 
* Hopkinson (2), 147. 
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principle of reflection, be analytically continued to every point of 
the lower half-plane. The function is thus seen to be regular at every 
point, including infinity, and so is a constant: its value is equal to 
that at infinity, viz. zero. Hence 

a i v(l—aj) | V(l— aj) V(1—u?) J —a,”) 
du J(1—u®)| u—wy u—t, u—U; u—t, | 





(3.3) 
AK 


C2kIl 
v-plane 
Yo 
cc 
Fia. 5 





It will be found simpler to transform the u-plane into a v-plane 


by means of the relation 
Qn 
20 


“= ——., 

1+v? 

The area of motion is then represented conformally on the upper 
half of the unit circle in the v-plane (Fig. 5). We find 


(3.4) 


Q = los | > See sate | Hose, (3.5) 


1 mt — , 
1—v,v 1—6,v v—v, v—3, 


the constant of integration being determined by reference to point 
B or D. Hence 

dw 1 (1—v, v)(1—#,v) (v—v,)(v—4) 

dz c (v—v,)(v—#,) (1—v, v)(1—d,v) 

Also, if vg denote the position of C in the v-plane, (2.3) gives 
dw B(L+v3) Lv (X—vy)(1—vj »)(v—8)(1— 6} 0) 
dv e(v—v9)(1—vy v) 1+ v? (v—v,)(1—2, v)(v—6,)(1—4, v) 

32 2 

oii b = ¢ tet e) 

(1+-7,°)(1+-@,") 


Define new constants, all real, as follows: 


(3.6) 


» (3.7) 


p = ||’, q = —R{r,}, r= Sry}, 
p= |y4/, g’ = —R{v}}. 
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Then dw (1+2qv+pv*)(v?+ 2q'v+p’) 
dz c(v®+2qv+-p)(1-+2q'0+p'r?)’ 


while, from (3.7) and (3.10), 
dz + b(1+-2) 1—v?2 (1+-2q’v+-p’v*)? 


dv (v—v»)(1—vp v) 1+? (1+ 2qv+- pr?) ’ 
1—p)?+49? 
and from (3.8) b=<¢ ( py + ae 
(1—p')?+-4q° 

The constants p, g, p’, g’ are not all independent, but are con- 
nected through (2.3). The relations between them can best be 


expressed by the identity 


(v2-+ 2qv-+-p)(1-+ 2qv-+po®) + ——", (v—v9)(1—v9 v)(1-+02) 


ot 
1+ 


l 
= —(v?+ 2q'v-+-p’)(1+2q’v-+p'v®), (3.13) 


where h = xr(1—p). (3.14) 
In the following investigation, the four basic arbitrary constants 
determining the motion may most conveniently be taken to be p, gq, 
Vp, k. The values of these being assigned, uw, follows, and it is neces- 
sary next to find p’ and q’. The identity (3.13) yields the following 
quadratic equation for p’: 


1+p* 


/ 


p 
(1—p)?+4q? , [{(1+-p+2q)?+4h(1—a)}{(1+-p—2q)?—4h(1+-u)}], 
2(p—huy) = 2(p—huy) 
(3.15) 
The root which is numerically less than unity must be taken (see 
Fig. 5). Then , pp g(itp)+h 
= l+p’ p—hu, — 
When v = i, = =: Le where 6 denotes the ultimate inclination 
dz c 
of the stream to Ox. Hence 
(1—p—2ig)(1—p'+2iq’) gas _ 2A —p)q’—(1—p')a] 


et —- . 7 > 7\? i aoe 
(l—p-+ 2%q)(1—p’— 219’) (1—p)(1—p’)+49q 
(3.17) 


(3.16) 
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Now use of the identity (3.13) in (3.11) gives 
dz __e(1+v8)(1—v®)(1+ 2g’ p'e%)(0® + 2v+-p) 
dv ~ (v—v )(1—v 9 v)(1+-v?)(1+ 2qv+ pe?)(v?+ 2q’v+-p’) 

2he(1—v?)(1+ 2q’v+p’v?) 

2\2/,.2 , *\° (3.18) 

(1+-2qv-+ pe*)*(v? + 2¢'u-+p’) 
This can be put into the partial fractions 


dz A, Agu 44s gut+A, 
1+v? 


dv v—W% 1—vpv 
2he| C,v+C, ra A,v+A, 4 (1—v?)(1+ 29’v+p’v?) 
1+2qv+pv? | v?+2q’v+p’ © (1+ 2qv+ pv?)*(v®+ 29’v+-p’) |’ 
(3.19) 








— 


$4. 


| 





where all the new constants are real, and 


Dey'95 4,,.2\2 2 , ‘\2 
hgon fer) hain bt eeate . (8.20) 
1+ 2qvo+ pr vt 2qvotp 


; 1—p’+ 2%q’\? 
—A 4 = 
on fue) 





1+ 2q'0,+p'vj 


<_(142mtmt lan 


C, 
oe ee 
mt F b(1—p) 
determining C, and C,. 

Making use of the fact that dz/dv has no singularity at either zero 
of v?+-29q’v+p’, we find that (3.19) reduces to 

dz A, Aygt A,v+A, 

dv ar pt aa Cai 1+v? rt 

1+-2qv+pr* (1+ 2qv+-pv?)? 


__ (vg—1)(v7+ 2q'r,+-p’) 
v2(1+ 2q'v,+p'v?) 








iy 


where f+0 s= 


and 


, , C, 
(1+29'e,+7'e)| (24 C,)op-+ar)+1C0 | 
| 


= (ef —1)(p'+-9'0y)— (07+ 29’) +P’)— V1 (1+-9'%)(C5+CG %4), 
determining C,, C,, C;, Cg. 


3695.10 
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Defining further real constants by the relations 


C.q—C; 


2r? 


C.p—C; q 


C, = - , G= 


¢ > 
2r? 


: CC; 


] 
C= F[GtGtO—G+G)t], Cy = FE, (8.28) 


we obtain 


dz A, = Agu% Asvt+ Ay | 
dv v—v, 1—wv 1+? 
4 hel GEGMetalp)+2rCy , d (_CyvtQ, \) 
1+ 2qv-+ pr? dv\1-+-2qv+ pv? 
Integrating this, we have 
z-+ constant 


= A, log(v—v,)+ A, log(1—vy v)+ $A, log(1+v?)+ A, tan—v+ 


+ 2he 


1 ty i 
coe $+2¢ytan-P7 54 + Clog(1-+-290-+po%) | 

The functions tan~!v, tan-{(pv-+-q)/r} are uniform within the semi- 
circular area (Fig. 5). Take the branches, the values of which at the 
origin are, respectively, 0 and (pr.argv,— 47). The functions may 
then be written in the alternative forms 


a 4 7—i(pu+ os 
3i log 1 ° and lilog’ Upotg) (principal values) 
Tw 


-r-+u(pv+q) 
respectively. Take the principal branches of the other logarithms 
also. 

When v = +1, z = +(d+4iz). Subtracting the equations result- 
ing from the substitution of these values, we find A, = 1. Hence 


p= (Lbemotons 


ate, lana 
1+ 29¢'v9+p'r} 


ae (v§+ 29'v9+-p’)(1+ 2qv9+ pro) and A. 


> 


52 ‘ 9 | Day's 42 
(ep + 2Qvo +P) 1+ 2q'v9+ pvp 
Also, by making use of a relation derivable from the identity (3.13), 


—A,+iA, = 2ce’®, (3.25) 
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Thus the relation between z and v is 
zta+thin 
= log(v—vy)+c*log(1—v, v)—c cos 8 log +e 


; + 2c sind tan-lv + 


Te amet? 2 
+2he| oer Tt +20, tan ~ + C,) log(1+-2qv+pv*) |. 

(3.26) 
where, by adding the equations for v = +1, 


_ 14+ 
ae 


C,(1+-p)—C;, 2q¢ _pt+4q 1P—-4 
+2he| (Ip) 4g? +0,}tan i —tan pan 


a log(1—v2)+ 





+1Cylog{(1+p)*—445}, (3.27) 


and, by subtracting the equations, 


2 } 
— log St %_ dnc sin 6 + 
2 1—v 
C, 2q—C(1+p) 
(1p) —4¢? 


d= 





+2he| —6,{tan~ Prd + tan-P—a re 
rT 
1+p—2q 


4C,, lo, 
Te Tp +a 


| (3.28) 


4. The free stream-lines and their asymptotes 

The free stream-lines in the z-plane correspond to the curved 
boundary of the semicircle in the v-plane (Fig. 5), on which v = e®, 
the lower one corresponding to the part on which 0 < @ < }m and 
the upper one to the part on which $7 < 6 < =z. 

Defining two real functions f,(@), f.(@), continuous on the semi- 
circular are, by means of the relation 
e+ 4 


f,(0)+tf2(0) = log(e”—v,)+c?log(1—v, e"®) —c cos § log : 


a 


ae reg C,e%+C 
" id 9 ee 3 
+ic sin 8 log(e +i)+ 2he| eat 


r—i(pe+q) 


So 
+ #0, log r+i(pe”+-q) 


+ Colog( 1+ 24+ pet) | (4.1) 
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we deduce the equations of the lower and upper free stream-lines 
in the parametric form 


x+a-+c cos 8 log{2 sin(}a7—46)}—c¢(3a-+ $0)sind = i (A), 
y+4n—c(4a—}0)cos 8+ csin 6 log{2 sin(}7—}6)} = f,(0 


and 
x+a-+ccos 68 log{2 sin(46—4}7)}—c¢(30— i 


y+4n+c(40—$)cos +c sin 8 log{2 sin(}@—47)} = 
respectively. The equations of the lower and upper asymptotes are 
respectively 
{x-+a—f,(47)}sec5—ac tand = {y+ }4a—f,(47)}cosec 6, 


{x-+-a—f,(47)}sec 8-+-7e cots = {y+ 37—f,($7)}cosecd, (4.3) 


where 


fda) + fa(3 


= log(i—v,)-+c* log(1—v, 7) — $act cos 6— $e sin 8+-ci sin 8 log 2+ 


+2he| Cri+C, 


1—p+2iq +? iC, log— Pa Cy log(—p+2ia) | 


Putting the constant in equation (2.2) equal to —}im we have 
w = log i(uy—u) —ix|log(|u—u,|/|w—a,|) +7 arg(u—u,)—i arg(w—1%,) |. 
(4.4) 


The second argument is one-valued in the area of motion and so 
can always be given its principal value. The first, however, is many- 
valued. 

On the free stream-lines, v = e’® and u = sec@, so that 


w = log i(uy—sec @)+-«[arg(sec 0—u,)—arg(sec 0—%,)]. 


The lower free stream-line corresponds to 3{w} = —4m and the 
upper to 3{w} = 47. On both 


R{w} = log|sec 0—u,y|+-«[arg(sec 0—u,)—arg(sec 0—%,)]. 


If s denote distance along the stream-line measured from D or B, 


s = clog es +x ang(“Soy ane a ")]. (4.5) 
—Uy —Uy 


+1 +1 +1—4, 
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Also, if % denote the angle of contingence at a point of the stream- 


line, then d 
D ap te 

ce'¥ = ; 

() if 


whence et — (29+ 2ge"?+-p)(1+-24 e+ p' *) (4.6) 
(1+ 2ge'® + pe9)(e219+ 29'e! +p’) 
The intrinsic equations of the free stream-lines are thus obtained 
in terms of the parameter 0. 


5. Determination of constants associated with the vortex 

The position of the vortex being denoted by z,, we have, putting 
z= 2%, v= v, = —q-+7r in (3.26), 
2+a+hin = 

log(—vy—q+tr)+c? log(1+-v) q—tv,r)—c cos 8 log 





1+q?—r?—2igr 
+9 ; ar 





Ltrtt sel —Crg Og iGy 
0q (1—p)(1+-p—2q?+ 2igr) 


r(1+p)—tq(1—p) , 
+10, log — “isin =P} + Cylog{(1—p)(-+-p—24*+-2iar)} |. 


+c sin 6 log + 


Hence 
x,+a = 
} log(vz+ 2v) q+-p)+ $c? log(1+ 2v9 ¢+- v3 p) — $c cos 8 log 


1+p—2q?)—C(1—p)q _ 
1—p a —p)?+4r} 


—C, arg{2r?— p(1—p)—2igr}+-C,, log| (1—p),/{(1—p)?+ ar], 


“i  . 





—esin8 arg(1—p+2ig)+ 2he| a( 


and 
yt do = arg(—vy—q-+ir)+e% arg(1+-r9q—ir9r)— 
ae 
+ 4g} * 
Cil+p)—G.20r 6 ipl one 
(1—p){(1—p)?-+.4r9} © * —p 
ey 
(5.1) 





—c cos 6 arg(1-+-p—2r?— 2igr)+e mtg 


+2he| 


Principal arguments must be taken in all cases. 
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Now let U,, V, denote the components of the velocity of the vortex, 
if free, or of the fluid at the vortex, if bound. If we imagine the 
(bound) vortex replaced by a small two-dimensional body, e.g. an 
aerofoil, round which the circulation is 2«7z, there will be a force on 
the body, say of magnitude F, at angle 8, to Ox,, given by the Kutta- 
Joukowski law 


F, ef: — —i2knQ, e'% = 2nn(V,—iU,) 


(taking the unit of density to be that of the fluid). Further, 


oe ee dw\* , 
Fie B =3/ (=) dz 
(2,) 
=5 [ : 1+ 1—v? /v?+-2q'v+p’ ihe 
= J c2 (v—V9)(1—%% v) 1+ ye yet 2qv+p 
(v,) 


[ b(1+-vG)v* beers 1 


(1—v9 v)(v—v%) v?-+1\ 1+ 2qv+ pr? v 





(1)v,) 


dz 2 h “7y 1 FY 
= dv = — = (iCy+ Cr). 


2h 


Hence = = (Gro + 1 9) 


(Cy+ tC jo). 


and U,tW, = Qe = _ 2 
We shall determine f, so that 0 < 8, < 7, and allow F, to take both 
positive and negative values. 

Again, if w,, denote the potential due to the vortex in infinite 
fluid, and w, the complementary potential w—w,,, the radius of 
curvature FR, of the stream at the vortex is given by 


] 


23 _ 1 {priate 


hk, 2° \ dz? |, 
We find that 
w, = log(u—uy) +i log! 
2, 


Cn 
n! 


d"z 
where ‘= (n = 
7 
du"), 


(u—u,)” _ +ix log(u—4a,)—}iz, 
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This leads to the result 


_1sf, f a 2 —P\A+v7)*(Ag-+v, Ap+ vj gAp— vi pAg) 
1 be,(1+-2q'v,+p'vj)* 


tales) +a -mleta—a)|> © 


where Ap = p’—p and Aq = q'—4. 

To apply these formulae in practice, it is necessary to deduce the 
constants p, qg, Vy from the physical data d, x,, y,, x. The equations 
are tractable in three main cases: 





(i) if the vortex lies appreciably downstream; 
(ii) if the vortex lies near an edge of the orifice; 
(iii) if the vortex is weak. 
The solutions for these cases (with overlap zero) will be given later. 


6. Problem of stream with free boundaries only 
The solution for the case where the fluid is limited by free boun- 
daries only will now be obtained for general values of the vortex 
strength. 
A 





0 
2K TI 


a 











ra 
Fic. 6 

The ultimate width of the stream is equal to the initial width, 
7 units, say. Take the velocity along the free stream-lines as unit 
velocity. Let the strength of the vortex be 2«z. Let the angle 
through which the stream is turned be 56. 

Take the mid-line of the undisturbed stream as z-axis. Let the 
y-axis be chosen to pass through the vortex, whose coordinates are 
(0, y,). Let w represent the complex potential. 

We seek a conformal transformation « = u(z) which shall repre- 
sent the area of motion on the upper wu half-plane so that C,C’ 
transform to the origin, A, A’ to the point at infinity, while |u,| = 1. 
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Motion in the u-plane represented by the complex potential w is 
due to a vortex at u, of strength 2x7 and a source at the origin 
discharging fluid at rate 7, while the real axis is a bounding stream- 
line. Taking account of the image system, we obtain 

w = log u—ix log(u—u,)+7« log(u—4a@,)+ constant (6.1) 


dw | ix(ujy—%,) — (u—u;)(u—a;) 


and —_— = = . 
dus us (u—u,)(u—%,) u(u—u,)(u—4,) 
where uw; represents the stagnation point 

We deduce that 

U+G, = Uy +4,+%K(u,—a,) and |u,;|=1. (6.2) 

Now Q must be purely imaginary on the real axis of wu, be finite 
at infinity, and must have logarithmic singularities at w,, u;, with 
coefficients of the logarithmic parts --1 respectively. By an argu- 
ment similar to that used to establish (3.3) we find 


() 
dQ _ ] 1 ae 1 1 (6.3) 


—— as a ee = i+ =~/* 
duos u—-& Uh, Uy, U4, 


Integrating, we obtain 
dz an gis (U—%y)(u—%,) 
dw (u—a,)(u—u,) 
Putting wu = 0 in this we find 
= eis 
” = 9 
ti 
so that we may write 
. = eifA-28) uy el A+t8), 
where A is real. Equation (6.2) now gives 
2 sin Asin $5 
sin(A— 458) 


From (6.1) and (6.4) it follows that 
dz e®/u—a;\? 
du u \u—ai,) © 
Integrating this equation, we obtain the required relation between 
z and uw, namely 


K= 


z+a,+ia, = logu+(e®—1)log(u—a,)—(e"4?—1)2a, /(u—a,), 
(6.7) 
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where the logarithms take their principal values, and a,,a, are real 


constants given by 


a, = —2sin® {6 log{2sin(A—35)}—Jorsind + - int a 


a, = $n—2sin? 43{7—A+}45-+ tan }9}. (6.8) 


= 


We deduce that 
Y, = —(a—A+ 48)cos +sin 8 log{2 sin(A—}8)} — 


sin(?5—A), 


2 sin? 35 cos(A+-36) 
~ sin(A—48) 
(6.9) 
To determine the quantities U,, V,, F,, B,, we have, as in §5, 
F,e#: = —2ixnrQ, e'™, 


and F, eB: — _ al (a) dz. 
dz 


(2) 
These formulae lead to the results 


F, ef: = —2rsin 35 e(7 +902 
sin(A— }5)cos 33 eidle 


and U4, = Q,e = 
it 1 Q; sinA 

so that the direction of motion is equally inclined to the initial and 

final directions of the stream. 


We again have the formula 


1 
ae —75(e erin 
R, Q: 


and, after some reduction, obtain 


1 «sin*(A—}8) P P 
BR, = 3sin®Acos 8  cosec A—} cot?(A—}8)+ — 


fen 
d*w,\ 
dz* I 


3x cosA 


4sin?) 1; 


16 cant =i 
(6.12) 


From (6.1) and (6.4) it is seen that the equations of the free 
stream-lines are 


8 = log|u|+« arg(w—u,)—« arg(u—4%,), 
ita) — (U—%y)(U— Hy) (6.13) 
(w—t,)(w—u) 
where w is a real parameter, positive values of w giving the lower 


stream-line, and negative values the upper. 
It is easily shown that one of the free stream-lines is inflected 
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unless A = 47. The inflected line is the lower or upper according as 
A < }ror > fr. 

When the vortex is weak, i.e. « is small, we have the first approxi- 
mations caer er. ox Setn., ee ee 

a, = —x«7, a, = in, 
= U,—1 = }xtany,, AV == k, 
F, = —2«z, B, = n+, 
x = «(4+} tan’y,). 
1 


The equations for the free stream-lines reduce to 


sin y, +e® ee 
b= ei, (6.15) 
sin y,-+-cosh s 


7. Conclusion 

In the above work the units have been chosen so as to give certain 
parameters convenient values. The results can, however, easily be 
generalized. If the distance between the fixed boundaries (or initial 
width of stream) be k, the general velocity Uj, and the fluid density 
p, we have merely to put 7x,/k for 2,, 7y,/k for y,, «7/U,)k for x, 


U/U, for U, V/U, for V, 7F,/pU2 k for F,, and so on. 

The writer, in conclusion, wishes to express his indebtedness to 
Dr. R. C. Knight for supervising the above investigation and for 
valuable suggestions. 
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By N. SIMMONS (London) 


[Received 16 February 1939] 

1. Introduction 

IN a previous paper (4) I have discussed the problem of the influence 
of free boundaries on the field of flow around a vortex in two dimen- 
sions. The solution has an important application in aerodynamic 
work, where it is usual to replace an aerofoil by its equivalent 
vortex system, whether in two or three dimensions, and by investi- 
gating the motion of the fluid in its neighbourhood, to determine the 
necessary corrections to the lift and drag coefficients, due to the 
aerofoil being placed in a finite stream instead of in infinite fluid. 
We limit ourselves here to two-dimensional motion, the aerofoil 
being replaceable by.a single rectilinear vortex. 

The problem we investigate is concerned with an aerofoil or vortex 
in a stream of finite breadth issuing from between two parallel walls. 
In particular, as indicated in § 5 of (4) we can deal with three cases: 

(a) (i) if the vortex lies appreciably downstream; 

(ii) if the vortex lies near an edge of the orifice; 

(iii) if the vortex is weak. 
There is also the problem (b) of the aerofoil in the free stream, which 
has been shown to be soluble without restriction.* This is, of course, 
the limiting case of (a) (i) and I shall not discuss it in detail. The 
approximations in each case entail some considerable reduction, 
which, although different in the three, follows the same general lines. 
I shall therefore be content to deal in some detail with (a) (i) in § 2, 
merely indicating in $$ 3, 4 the type of approximation required for 
(a) (ii) and (a) (iii) respectively. The results for all three cases are 
generalized in § 5. 

We have observed} that Prandtl’s method is open to objection in 
its use of image systems. Also the solution thus obtained and 
discussed by Toussaint}, applying primarily to problem (6) is used 
for the more practical problem (a), the assumption being implied 
that the fluid emerges from the walls in a direction parallel to them. 
It is reasonable to assume this if the distance of the aerofoil from 


* (4) §6. t (4) §1. t (2) 302-5. 
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the mouth of the tunnel is sufficiently large. In fact, the effect of 
the mouth is found to be negligible only if this distance exceeds 
1-5 times the distance between the walls. Another assumption made 
in Toussaint’s discussion we leave to § 7. 

From the generalized results of §5, we deduce formulae for the 
corrections to lift and drag coefficients of a two-dimensional aerofoil 
(a) in the mouth of a wind tunnel, and (b) in a free stream, without 
these assumptions. We compare them with the formulae given by 
Toussaint* for the symmetrical case, and suggest that his results 
require modification. 

The notation employed is, throughout, that of (4). 


2. When the vortex lies appreciably downstream 
A first approximation will be obtained, but it will be necessary 
initially to find some constants to second approximations. 
In the notation of (4), |v,—¢| is small. 
Put v, = i(1—ee*). 
Equations (3.15), (3.16) of (4) lead to 


v; = i(1l—ee-™), 


where cos w’ = [(1—«?)cos? w+ 2« cos w sin w}}, 


sin w’ = KcOSw—sin w. 

We have, to first approximations, 

1—p = 2€ cosa, q = —esinw, r=], 

I—p’ = 2€cosw’, q = esinw’, h = 2ke cos w. 
Equation (3.17) of (4) yields 

§ = 2(w+w’). 

Also, from (3.24) of (4), 6=c¢= ]., 
We next evaluate the constants C,,, and deduce the following values 
of the physical constants: 

k = 2sin }(w+w’)cos }(w—w’)sec w, 

a = log(1—v?)—2sin?(w+w’)log 2+ O(e), 


l+v 
0 +0(6), 
1—vy 


* (2) 304-5. 


d = log 
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x, +a = log(1+v?2)—log «+w sin(2w+ 2w’)+ 


+2sin? }(w+o’ Sean + O(e log «), 
Y, = —w cos(2w+ 2w’)+sin(2w+ 2w’ )log(2 cos w)+- 
+2sin? }(w+w == +O(€), 
U4, = ees 3(w+w’) ) ehid, 
cos }(w—w’) 
Fi ths ae then ET gn, (2.7) 
cos }(w—w’) 
(Note that loge is not necessarily large.) 

For non-overlapping ends, d = 0, and so v) = O(e). Thus, sub- 
stituting w+w’ = 46 and putting w—w’ = 7— 2A, we obtain the first 
approximations 

__ 2sinAsin $5 
= seals ‘ 


= log? +(Jn—A-+-}8)ain 8-+-2in n2poSnte) 


n(—45) +2 sin? 45 log 2, 
= —(4n—A+ }5)cos 5+ sin 8 log{2 sin(A—}5)}— 
. A+48) 
_9sin? 5 CO8( 
sin* 3 n(A— 98)’ 
sin(A—}8)cos 48 
sinA 
oe sin(A—4}5)cos 45 
: sinA 
F, = —2«7rQ, = —2zsin 35, B, = 4n4+BB. (2.8) 

Comparison with § 6 of (4) shows that these results are the same 
as those for a jet with only free boundaries (when 2,, however, has 
no meaning). 

In practice x, y,, and « are the given quantities. Equation (2.8) (i) 
then gives an expression for A in terms of tand, so that, on sub- 
stituting for A in (2.8) (iii) an equation for 5 is obtained. This must 
be solved approximately: the other constants will then follow. 

For the intrinsic equations of the free stream-lines, (4.5) of (4) 
gives the first approximation 


= U+i% = 





a, = 38, 





8 = log|sec 4|, 
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so that, on the lower and upper stream-lines, 
e® = +sec, (2.9) 
respectively. Also from (4.6) of (4) 
iy — [cos @—e sin w—e cos w e~|[ cos 0-+-€ sin w’ —e cos w’e'9) 
[cos é—esin w—e cosw e9]| cos 0+-e sin w’ —e cos we? : 
Near the orifice, where «sec @ is small, this equation reduces to 
% = —4etan6@cosAsin }d. (2.10) 
From (2.9) and (2.10) we deduce the intrinsic equations of the lower 
and upper free stream-lines as 
yb = F 4e(e*— 1)! cos Asin 45, (2.11) 
respectively. In the neighbourhood of the orifice the jet is sym- 
metrical. 
For the asymptotes, we find 
fi(47) = —3asind, fo(47) = (1—cos 8)$2-+sin 8 log 2, 


so that the equations are 


(lower asymptote) y = #tand—4msecd+(2sin5—tan d)log 2 
log + 


(upper asymptote) y = xtand+ }7secd+(2sin5—tan6$) , 
(2.12) 


When the vortex is weak we have the first approximations 


l 
2, = log yg, = A—is = 46, 


=] ? < 


= a xX = K, 
F, = —2kxz =r —— «(4+ }tan’y,). (2.13) 
1 , a7; R, 314 1 . 
For an accuracy of 1 per cent., x, > log 100 = 37 approximately, 
i.e. the vortex must be situated at a distance from the orifice at 
least 1-5 times the width of the jet. 


3. When the vortex lies near one edge of the orifice 
The vortex is now near to B or D. Choosing the former without 
loss of generality, we have the position v, such that |v,+1{ is small. 
We may put pee = 5-t-a0¥, (3.1) 
with the resulting values 
p = 1—2e cos €+€?, = 1—ecosé, 


r= ening, = xe*sin 2€. 
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We find that « = O(c?) and then proceed in much the same way 
as for the preceding case, retaining where necessary higher approxi- 
mations than the first. Omitting the rather laborious reduction, we 
find that, if we define 


a= 1—see¢| }008 26 +5 (1—Ssinzg)? 
B = tané, : (3.3) 


the position, velocity, etc. of the vortex in the z-plane may be 
represented by 








my 2p i ar 
x a@—a app py TAT aT 
1 
— 3af(1—«)?— B? 
+a | Barf (1 —or)?—B7}+ 


+ $a(1—2a-+ 0? + B?— 208?) — 4a8? log ro 











—2{a8—2a%(1— BP) +-a(1-+ 8) — 264 1+ pr) A, 
= ®,(a, B), 
37—Y, 28?(4a—1) 
x  af®—a){(I—a) +BY 
1 
af{(1—«)*—?}+ 8a tan—1p— 
+a-a aL {( )?—B?}+ 808 tan—1p 


— 2a? —208(1—f8)+a( 1+) — 21+ pr) PEEP) 


=: ®,(a, B), 
a 1—2a-+02— 462+ 6oh?—a262— pt ~ O(a, 8), 


T ati = 
kes (a—2){(1—a) +B} 


, 4a 
R= oe = Mab 
F, cos B, = 2«7®,(a, 8), 
F, sin 8B, = —2«7®,(q, ). (3.4) 
Also 8 = [2«®,(a, B)|!, 
and the equations of the lower and upper free stream-lines at finite 
distances from the points of departure are 
% = 8(tanh 4s)}, x = 8(coth 3s) (3.5) 











respectively. 
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In any practical case x,, }a—y,, « are the given quantities: «, B 
may be treated as parameters, and a table of corresponding values 
of x,, Jn—y,, U,, V,, ete., constructed; 8 must be positive, while « is 
cosignal with x. Also the inequality « < 1 is useful when « is positive. 


4. When the vortex is weak 

In this case « is small. The process of approximation requires 
the constants to (at most) the second order, and these are readily 
obtained. We take v, = O(x) and find that this includes the case 
where the boundaries are non-overlapping. The values of the physi- 


cal quantities for this case, deduced from §§ 3-5 of (4) are found to be 
AU = U,-1 
« | sinh x. ts a aia o+sinh x,\4 
= 5-| sinh2, tany,+ —siny, cos y,—2 sin y,|——___ 
2c o o-+-cosh x, 
: « cosh x sal 
AV = = ———(e++sinhz,), — (4.1) 


Iq? 


where 
= |coshz,| = (cosh*x,—sin*y,)! = (sinh?x,+cos*y,)!, 


kcoshz,, |. 
= —2kn, B, = in-++——— (c++ sinh 2), 


€ 2 
207 





3 8 


| , 3 ( __ 2sinh x, sec*y, 4 sinh*x, cosh, —sin?y, —H)| 
4 3 


CG Co 
K If — “ ! a 
$= <Vi(o+sinh x,)(o-+ cosh x,)}. 


The equations of the lower and upper free stream-lines become 
respectively 

_ $6(1—e-**)? e*-*1(a-+- cosh x,)+e” sin y, 

~ o+cosha, cosh(z—2,)+siny, — 





In the symmetrical case, i.e. when the vortex lies on the mid-line, 
the equations become 
AU = 0, AV = 3x(1+tanh2z,), (4.2) 


> 


By = 407+4«(1+tanh2z,), 


R, = «[4—}(1—tanhz,)?], 
8(1—e-**)4 


= x,/{2(1+tanh 2,)}, y= pena” 


i= 
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Here, 5, 4, AV, R,, A, and 8,—4z must be odd functions of «x, and 
so the corresponding results are actually true to a second approxima- 


tion. 


5. Summary of results with full dimensions 
In the preceding sections the units have been chosen so as to give 
certain parameters convenient values. If the distance between the 
fixed boundaries (or initial width of stream) be k, the general velocity 
U,, and the fluid density p, we may generalize the results by putting 
wa,/k for 2, 7y,/k for y,, x7/U,k for x, U/U, for U, V/U, for V, 
aF,/pU2k for F,, and so on. We obtain 
(i) (The vortex appreciably downstream). 
_ 2U,k sin Asin $5 
~~ -sin(A—38)’ 





: log + + +(}7—A+} jo)sin3+ 2sint}s ATH) + osint45 log 2 ; 
7 € n(A—48) 


a 
919 COS(A+ 48) | 


2 sin?} 


sin(A—}58) 


sin(A—}8)cos 45 
sinA 





Q, Uy 


> 


F, = —p2xrQ, = —p2rsin $5, 


When the vortex is also weak, 


aca 
AU 


a 


al ¢ T KT ~ 
Ff, —p2xrU), By = 3 U, og (5.2) 
The formulae, apart from that for x,, apply also to the case of a 
vortex in a stream with only free boundaries.* 


* Cf. (4) § 6. 


2695.10 x 








306 N. SIMMONS 


(ii) (The vortex near an edge of the orifice). 











The lower and upper free stream-lines are respectively 


»\1 a 
b= a{tanh 7) and = A(coth TT , 


2k 2k 
Also +a, = 77 91(a:8): +hk—-y, = 7 ®,(«, B), 
0 0 


according as the vortex is near the lower or upper edge. Further 
U, = U,%5(«, B), yi +U,®,(«, B), 
F, cos B, p2x7U,®,(«, 8), F, sin B, F p2x7U, (a, B) 


(the signs being determined as above). 


The functions ®, (r = 1, 2, 3, 4) are those defined by equations (3.4). 
(iii) (The vortex weak). 
~ | 


o = |cosh— 


> 


C ° am & 39 
5 = ral (e4 sinh E(0 i cosh) |* 














es aed 72x i. iia 
- inate em(x—r Iki gt cosh —) + e7k sin 72 
6(1—e-272/k)s k 


b= 2 








TX 7 ~ ae 
o+cosh —} cosh 7 (x —a,)-++sin Wi 


v 




















: KT | . , 72x 7 l . m 1 
AU = sinh - Itan “A sin “#1 cos 741 — 
2ok 7 iy oC WC V 
. . We 
(o+-sinh =*\ 3 
> ae . 
ein 34{ 3 
: Tx 
” \o+cosh—! 
+ v 
y KT Ps TX TX 
AV = ——|o+sinh — |cosh —, 
20° 4 - 
a SIRE 9 TY 
fm, 4 2sinh —~ sec? =} 
KT7IT* ‘y 
—_— = “ ome =" 
r i 
Rk, U, k?| 3 ( o 
. ao AL Te . ot ey) 
sinh? —! cosh? —! _gin? 771 cos? 71 





| R k k k \| 
as o ’ 
R= 


—p2«rU), B, = 30+ = (5.4) 
0 
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When the vortex is symmetrically placed, the formulae are 


__ p—27z|k\4 
. - Ty (2+2tanh $= ai --¢ ee 








Uk 1p etmei—aik” 
KT TX 
gee Fs ae 1 | —- 
AU 0, Al 4 +tanh E |, 
1 km? fl 1 ara,\? 
Le a= ie 
R, rl 5 3( a” i} 
Fm peat, Ds tack or (5.5) 
0 


All these equations, except that for AU, are correct to a second 
approximation. 


6. Application to aerofoil theory 

There are two applications: (a) to an aerofoil in the stream of a 
wind tunnel, (6) to an aerofoil in a free stream of limited breadth 
extending to infinity both upstream and downstream. 

Putting I’ for the circulation in the clockwise sense, we have 


I’ = —2cr. (6.1) 


In all practical cases xx/kU, is a small quantity, so that the problem 
falls within the purview of case (iii) (5.4), or of the case with free 
boundary and weak vortex.* 

The stream is deflected downwards, say at an angle y (= —8). 
Let L’ denote the lift when the inclination of the direction of zero 
lift to the walls in (a), or the initial direction of the stream in (bd), 
is x, and let C;, denote the lift coefficient of the aerofoil corresponding 
to an angle of incidence a. Applying the Kutta-Joukowski theorem, 
we have L’ = pI’ U,, or, to a first approximation, 


L’ = pI’U,, (6.2) 
This relation is not true exactly, as is the case in unlimited flow. 


Now the effects of the limitation of motion on the flow around 
the aerofoil are 

(i) the angle of incidence is increased by A,«, say (= AV/U,): 
A, « is of the same order of magnitude as a, and so the correction 
due to this source may be considerable; 


* (4) §6. 
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(ii) the velocity of incidence is increased by amount AU: in the 
symmetrical case AU = 0; even in the asymmetrical case AU/U, is 
of order x, so that the effect here is negligible compared with that 
in (i); 

(iii) the stream incident on the aerofoil becomes curved. The 
curvature gives rise, according to Carafoli (1), to an increase in 
circulation of amount equivalent to a virtual increase in a given by 


A,a = (6.3) 


c 
4R,’ 
where c is the chord of the wing. 
In practice, A, « will be considerably greater than A,a, but the 
second may not be negligible. We shall write 
Aa = —A,a—Aga, (6.4) 
where Aw will now be a positive decrement. It is known that, to 
a sufficient order of accuracy, 
C,/« = 2s, . (6.5) 
a constant for the aerofoil. Hence the lift coefficient to be taken 
into account here is 28(a— Aa), and 


L’ = }pe2s(a—Aa)U2 = jpet,(1—"*) U2. (6.6) 
a 


Defining an apparent lift coefficient C;, so that 


L’ = 3pc0;, U3, 


and 
we have 


or 


giving the true lift coefficient in terms of the observed quantities. 
We also find, from (6.2), that, to a first approximation, 
, , T 2 
F — 40,cU). (6.10) 
Also, if AC, denote the correction to the observed drag coefficient 
taken as a positive decrement, 


ACy = Cy, cot Bj = —C_A,a. (6.11) 
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We proceed to apply these formulae: 
(a) Aerofoil in wind tunnel. 


r 
eS Cae ee. 


(Symmetrical case.) ag formulae of case (iii) give 


= i ( +5 Rai ny, 





== B= -=---9 





k 
1+tanh 
(stan) 


| zy 2 
A= 3(t-+tanh 7) 4 [1—§ (1—tanh =) | (6.12) 
With the usual arrangement it will generally suffice to take 


hon 3(1-+tanh 2) + (6.13) 


(Asymmetrical case.) The formulae of case (iii) give 
Cre 5 WE, ma,\ |4 
—- [(e-+sinn E }(o-+eosh E )| 2 ( ) 
A, a = _ Cre o-++sinh =~! cos sh ™=1, 
~~ Sat k &’ 


TY 
k 





+ 
o 


2sinh 1 sec? 71 
( , k 
sinh? cosh? 2! — sin? 1 cos? 7 


eet i) 
Co 


( - 2 sinh “sec? A 








k k 8 si 


oc 


53] ¢+sinh 7 [cosh cs | + 73 ial -3 


4 : k 





‘ TX TX . oF 
sinh? . cosh? _ —sin? A 1 oog2 71 | 








310 N. SIMMONS 


With the usual arrangement it will generally suffice to take 


i 31+ ta nt) +5 sin? 24(1—tanh +f (1+ jtant 7}, 


k 12k\ "4 
(6.15) 


(b) Aerofoil in a free stream. 





K-- Bp --------> 





7. Conclusion 


The results deduced by the above method for the case of an aero- 
foil placed symmetrically with respect to the wind tunnel are 


C, atx dite — a} 


Cr 4k k ° 6k 


and ACp = a (14-tann 2), (7.1) 


8k 
These results are accurate to the second order. The corresponding 
results from Prandtl’s method, given by Toussaint,* are seen to be 
of entirely different form. This appears to be due chiefly to the 
substitution, at one point in that analysis, of 3 for C;/2a. We suggest 
that s = C,/2a leads to a better result. 

It will be convenient to deal here also with another point, referred 
to in §1. We note that no allowance is made in Toussaint’s work 
for the effect of the curvature of the stream at the aerofoil, although 
the correction for this is found to have weight in the case of a wind 


* (2) 4.21. 
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tunnel of closed working section. The present analysis shows that 
A, «/A, « is of the order c/4k, so that, although the effect of curvature 
is small compared with that due to the alteration in the incident 
stream, it is not always negligible. 

To compare the numerical values of the corrections given by the 
two methods, we have the following table, calculated for the typical 
case where c/k = 0-20, and for various values of 2,/k. 


-~=010 020 030 040 050 0-60 oe) 


= 1-243 1-287 1-318 1-337 1-349 1-355 1-363 


Present 


method 
—? — 20-033 0-039 0-043 0-046 0-048 0-049 0-050 


Prandtl = 1-117 1-238 1-341 1-415 1-442 1-491 1-527 
and 
Toussaint ae = 0-017 0-034 0-049* 0-060 0-066 0-071 0-076 
L 

It is seen that the latter corrections are in deficit for small values 
of x,/k and in excess for larger values. The discrepancy ultimately 
increases with x,: it also increases with c/k. The point at which the 
formulae for C;, agree is approximately the same as that at which 
those for AC, agree. The formulae proposed above do not lead to 
zero corrections when x, = 0 as do those of (2). 

We note further, in our treatment of problem (a), that, provided 
that the unabridged formula for A be used, there is no restriction 
on the value of z,. It is thus possible to put x, = —0oo, when the 
stream is entirely enclosed between two fixed walls. Corrections are 
thus obtained for the wind tunnel of closed working section, which 
in the symmetrical case are in agreement with the results of the 
Prandtl method, amended by the substitution of s for C,/2« instead 
of C;/2«. We find the increase of lift quoted by Toussaint [| (2), 301] 
and not the decrease obtained by Sasaki (3), who, however, replaced 
the aerofoil by a flat plate. 

It is again a pleasure, in closing, to acknowledge the supervision 
of this work by Dr. R. C. Knight, and the advantage of his sug- 
gestions. 

* The numerical results for the drag, given by Toussaint [(2), 305], do not 
appear to follow from the formulae given. 
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ON THE LOBATCHEWSKIAN TRIGONOMETRY 
OF A STATIC SUBSTRATUM 
By G. J. WHITROW (Ozford) 
[Received 7 September 1939] 

1. Introduction 

Ir is generally agreed that the most novel and significant contribu- 
tion which Milne has made to cosmological theory is his discovery 
of the two families of time-scales, ¢ and r. 

The family of t-scales comprises all clocks, kept by any particular 
fundamental observer, in terms of which the substratum is described 
as uniformly expanding from a point-origin at the singular epoch 
t = 0. The family of 7-scales, on the other hand, comprises all clocks 
kept by the same observer, in terms of which the substratum is 
described as static. Given any one t-clock Milne has shown that we 
can associate with it a unique 7-clock which both reads the same 
epoch and beats at the same rate as the ¢-clock at an instant fp, 
which we can identify with the present. The relation between these 
two clocks is given by* 

T = ty log(t/ty) +t. (1) 

It was further shown by Milne that, if the observer using a parti- 
cular ¢t-clock plots his observations in a private Euclidean space and 
then transforms by (1) to a 7-clock, he generates a public hyperbolic 
space of constant negative curvature —(ct))~*. This result was 
proved by constructing the metric of this space. 

It is the object of the present note to obtain this result in a dif- 
ferent manner. I do this because I wish to point out an interesting 
parallelism between Milne’s discovery of the existence of a public 
position-space, for relatively stationary observers, and the apparently 
neglected theorem, due to Sommerfeld and, independently, to Vari- 
éak, establishing the existence of a public kinematic-space for ob- 
servers in uniform relative motion.{ Either result can be derived 
from the other by means of a simple formula, due to Milne, cor- 
relating a certain function of the speed assigned to a particle by an 

* Milne, Proc. Roy. Soc. A, 158 (1937), 324. 

+ Ibid. 159 (1937), 171. 


{ Sommerfeld, Verh. deutschen Phys. Ges. 11 (1909), 577; Variéak, Phys. 
Zeits. 11 (1910), 93, 287, 586. 
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observer using a given ¢-clock and the distance assigned to the same 
particle by the same observer at the corresponding instant according 
to the r-clock given by equation (1). For completeness I first establish 
this formula de novo and correlate Milne’s transformations with cer- 
tain transformations due to Lewis. This work is essentially one- 
dimensional because we need only consider an event in alinement 
with two observers. I then pass to the general case and give a new 
proof of Sommerfeld and Vari¢ak’s theorem, which is more concise 
than either of theirs because it is based on the use of 4-vectors. In 
this way I show how the distances and orientations of any event, 
according to two observers using congruent r-clocks, can be correlated 
and I also establish the invariance of the 7-epoch assigned to any 
event by these observers, without making the usual appeal to the 


geometry of ‘cross-country’ light paths. 


2. Two equivalent observers and a collinear event 

It is well known that the Lorentz transformation correlating the 
distances and epochs assigned by two observers,* using congruent 
t-clocks, to a collinear event can be regarded geometrically as the 
rotation through an imaginary angle of the axes of space-time co- 
ordinates. This is because the transformation can be written in the 
form 

’— x cosh w—ctsinhw, ct’ = ctcosh w—2x sinh w, 
1 Vie 


where cosh w = - sinh w = —_ 
(1—V?/c? 


(1—V*/c2)! 


and hence tanhw = 


The term w = tanh-!V/c has been called the relative rapidity of the 
observers and was introduced by Robb} because, in the Special 
Theory of Relativity, the rapidities, and not the velocities, of col- 
linear uniformly moving observers are additive. 

Recently, Lewist has drawn attention to an alternative way of 
regarding the Lorentz transformation as a rotation of axes. Lewis 


* In uniform relative motion of speed V, from coincidence at zero 
epoch. 

+ Robb, Optical Geometry of Motion (1911), 9; Geometry of Time and Space 
(1936), 406; Eddington, Mathematical Theory of Relativity (1930), 22. 

t Lewis, Phil. Mag. 20 (1935), 1096. His notation has been slightly altered 
so as to lead to no confusion between his symbols and mine. 
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introduces two auxiliary variables o, ¢, analogous to polar coordi- 
nates, and writes the transformation in the form 
z= cosinh d, t = ocosh¢d (5) 
x’ = cosinh(d—w), t’ = ocosh(¢d—w) 
where w is the relative rapidity. We see immediately that 
o = (t@—22/e2)}, (6) 
t+2/c 


lo : 
8 t—z/c 


x 1 
and = tanh-!— — — 
¢ ad 2 


The parameter o is invariant, but the parameter ¢ transforms into 
¢’, where i «= be. (7) 

Milne* has shown that, if (x,t) are the distance and epoch assigned 
to an event by an observer using a ¢-clock and if (A, 7) are the distance 
and epoch assigned to the same event by the same observer using 
the 7-clock given by (1), then the following relations exist: 


t = t,e7—lle cosh a 
Clo 


ain 
= ct, eto sinh 2 


Clo 


Consequently, Lewis’s parameters o, ¢ become 


c= (t2—2?/c?)t — é. eft —to)lto 
(9) 
x A 
¢@ = tanh-!— = — 
ct = chy 
Furthermore, since on the t-scale the observers’ relative speed of 
recession V is given by 
7 _ 2& A 
Vv =— = ctanh—, 
t Cty 
it follows that, if, on the 7r-scale, the same observers’ constant 
distance apart is A, then, as Milne has shown, 


v= ctanh. (10) 


Cho 
Consequently, if w is the relative rapidity of the observers using the 


t-scale, then A 
w= —., . (11) 
Cty 


* Milne, Proc. Roy. Soc. A, 159 (1937), 177. 
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We thus deduce from (7), (9), and (11) that, if the distances assigned 
to an event on their common radius vector by two observers, distant 
A apart, using congruent 7-scales are A and X’ respectively, then 


’ = d-A. (12) 


Consequently, the 7-scale distances between the observers and a col- 
linear event are additive. Similarly, from (9) it follows that the 
invariance of o implies that the respective conventional epochs 7 and 
r’, assigned by the observers to a collinear event, are identical. 


3. Two equivalent observers and any event 


We now pay attention to orientation and consider two observers 
and a non-collinear event. 

First, we observe that the invariance of the conventional epoch 
assigned by two observers, using congruent 7-scales, to a collinear 
event can be immediately extended to the case of a non-collinear 
event, without making any appeal to the geometry of ‘cross-country’ 
light paths.* This is possible because we can always locate a funda- 
mental observer at any event in the substratum. The conventional 
t-epoch assigned to the event by any other fundamental observer is 
identical with this particular observer’s own 7-scale reading at the 
epoch concerned. Consequently, all fundamental observers assign 
the same r-epoch to any given event.t 

We now consider the relations between the radial distances and 
orientations of any event relative to two observers using congruent 
z-clocks. Milne has found that there exists a public hyperbolic space 
of constant negative curvature —1/(ct,)*, it being assumed in the 
first place that each observer using his original t-scale correlates his 
observations with those of the other observer by means of the three- 
dimensional Lorentz formulae. We now see, in view of relation (11), 
that Milne’s result is formally identical with the theorem of Sommer- 

* For the necessity of such an appeal, or its equivalent, to establish the 
Lorentz formulae in space of two or more dimensions vide Milne, Relativity, 
Gravitation, and World-Structure (1935), § 36. 

+ Incidentally, if we transform back to the ¢t-scale, given by inverting (1), 
we find that the well-known invariance of the quadratic form, X = t?—P?/c?, 
where ¢ is the conventional epoch assigned to an event at P, according to 
any member of a pencil of equivalent observers uniformly expanding from 
a point origin at zero epoch, is independent of the validity of the three-dimen- 
sional Lorentz formulae and is purely a consequence of the one-dimensional 
formulae. This point is not usually made clear in writings on the subject. 
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feld and Vari¢ak, that in the Minkowski diagram the trigonometry 
of rapidities is Lobatchewskian. I proceed to establish this theorem 
de novo, by means of 4-vector analysis. I find it convenient to con- 
sider three observers, instead of simply two observers and an event. 
As we have already remarked, we can always locate a fundamental 
observer at any event in the substratum. 





Let A, B, C be any three non-collinear members of a pencil of 
equivalent observers in uniform relative motion from a point-origin 
at zero epoch, each observer plotting his observations in a Euclidean 
private space and correlating his observations with those of the other 
observers by means of the three-dimensional Lorentz formulae. Let 
the constant angle subtended at A, according to A, by B and C, be 
Q,; the constant angle subtended at B, according to B, by A and C, 
be Q,; and the constant angle subtended at C, according to C, by 
A and B, be Q,. Let the relative rapidities of A and B, B and C, 
C and A be wo, w;, w, respectively. Then the triangle, constructed 
in Fig. 1, with sides wy, w,, w, and angles Q), Q,, Q, will obey Loba- 
tchewskian trigonometry if 


cosh wy = cosh w, cosh w,—sinh w, sinh w, cosQ,, (13) 


sinQ,  sinQ, _ sinQ, (14) 


and if : = — = — : 
sinhw, sinhw, sinha, 


If V,, V,, V2 be the relative velocities of A and B, B and C, 
C and A respectively, then, by formulae (3), 


1 
(1—Vo/e?)” 


‘ V 
cosh wy = sinh wy = a Wea (15) 
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and similarly for cosh w,, sinh w, and cosh w., sinh w,. Consequently, 
1—V,.V,/c? 
(1— V2/c?)}(1— V2/c?)} 
(16) 
The expression on the right-hand side of this identity is simply the 
scalar product of the two four-vectors, 


( o/€ Hi) and fs Fab 
Yi ¥} esas 


where Y} = (1—V3/c?)!, Y} = (1—V?2/c?)! 


cosh w, cosh w,—sinh w, sinh w, cosQ, = 


It is therefore invariant under all Lorentz transformations. If we 
regard A and C as fixed observers, it follows that this scalar product 
has the same value howsoever we choose B from the pencil of equi- 
valent observers considered. In particular, if we identify B with C, 
this scalar product reduces to (1—V3/c?)-?. It follows that 
1—V,. V,/c* (17) 
= . j 
[2 [p2\s 2/p2\k? 
(1—V6/c*)*(1—Vj/c)! 
and hence that 
cosh w, = cosh wy cosh w,—sinh wy, sinh w, cosQ,. 
Formula (13) is, therefore, established. 
To verify (14) we observe that, since 


sinh wy 


it follows that sinh w,sinQ, 


(1 —V2 c?)t’ 


where v, is the component of V,, orthogonal to the direction AB. 


Since V,/c 1 ana (Vee 2 
ip a yi’ i 


are 4-vectors, we see that, if y, y’ are the coordinates assigned to C 

by A, B respectively in the direction orthogonal to AB, then, since 
= 2. 

* Y% U2 


18 
(1—V3/c?)t = (1—V3/c?)#?” (18) 


where v, is the component of V, orthogonal to AB. Consequently, 


sinh w,sinQ, = sinh w, sinQ,, 
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sinQ, 
and hence ape UE eee, 
sinh w, 


In this way formulae (14) can be established. 

We have thus proved that a substratum of equivalent observers 
uniformly expanding from a point source possesses a public kine- 
matic space, since the trigonometry of rapidity vectors is that of 
a hypersphere of radius 7. We can draw the following conclusions: 


(I) Any two rapidity vectors* can be compounded by the triangle 
rule of hyperbolic space. Incidentally we note that these rapidity 
vectors are not free vectors but line vectors. There is consequently 
no parallelogram rule. 

(II) The resultant of any number of rapidity vectors arranged in 
a chain in hyperbolic space is given by the geodesic drawn from the 
beginning to the end of the chain. 

(III) Finally, by equation (11), we see that, when the observers 
regraduate their congruent clocks from ¢-time to 7-time by formula 
(1), there exists a public position space, which is hyperbolic and of 
constant negative curvature —(ct,)-*. The relations between the 
radial distances and orientations assigned to any event by any two 
observers of the system are thus given by the appropriate formulae 
of Lobatchewskian trigonometry. 

* As we are primarily interested in equivalent observers, we have considered 
only rapidity vectors corresponding to motions from a common point origin. 
The proof, however, can be applied step by step to rapidity vectors in general 
as in the Special Theory of Relativity. 








ON A CONJECTURE OF GILLIS 
By 'T. VISAYARAGHAVAN (Dacca) 
[Received 8 November 1939] 


Suppose that f(x) is continuous and periodic and that the line 


y = mx-+c, where m and ¢ are finite constants and m ~ 0, meets 
the curve y = f(x) at a finite number of points only. Let C denote 
the curve y = f(x), and L the line y = ma-+c. I prove below that 
there exist finite constants m’ and c’ (m’ ~ 0) such that y = m’x-+c’ 
meets C at one point only. This result confirms a conjecture* of 
Gillis that, if C be such that every non-vertical straight line meets 
C in at least two points, then every such straight line meets C at 
an infinity of points. The existence of continuous periodic functions 
that give rise to such curves has been established by Gillis. 
Firstly, let m > 0. It is plain that LZ is on both sides of C, and 
ultimately is above C. Since by our hypothesis we have that LZ 
meets C at a finite number of points only, it will follow that there 
is a number €, and a positive number 7» such that 
fié) = mé+e, 
> m(é—5)+c = f(€)—més, 
+8)+e = f(€)+méd 
whenever 0 < 6 < 7». 
Now, let H be the difference between the upper and lower bounds 
of f(x). (4) 
We consider y = m’x-+-c’ where m’ and c’ are determined so as to 
satisfy the following conditions: 
m' >m > 0, 
m'n > 2, 
mé+e' = fié). 
We easily verify that 
m'(E—8)-+-e" < f(§—8) 
and 
m'(&+8)+e' > f(E+8) (9) 
whenever 5 > 0. From (7), (8), (9) it follows that the straight line 
y = m'x-+c’ meets C in one point only. We argue similarly if m < 0. 


* See above, p. 154. + Loc. cit. 151-4. 
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